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The  thermodynamic  pressure  of  superfluid  '*He,  for  seven 
densities  from  near  saturated  vapor  pressure  to  11  atm. ,  has 
been  measured  from  0.4  K  to  1  K,  using  a  capacitive  strain 
gauge.   A  comparison  of  these  measurements  with  calculations 
for  the  pressure,  based  on  the  phonon-roton  model  for  super- 
fluid  helium,  led  to  the  conclusion  that  anomalous  phonon 
dispersion  exists  in  helium  at  low  densities.   The  magnitude 
of  the  anomalous  dispersion  is  in  good  agreement  with  that 
deduced  from  specific  heat  measurements.   The  conjecture  of 
a  quadratic  term  in  the  expression  for  the  dispersion  curve 
does  not  agree  with  these  data.   Contributions  to  the  pres- 
sure not  accounted  for  by  phonons  and  rotons  were  observed. 
These  contributions,  probably  maxons,  increased  in  magnitude 
and  became  noticeable  at  a  lower  temperature  as  the  density 
was  increased. 
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CHAPTER  I 
INTRODUCTION 

Historical  Background 

Liquid  helium  has  been  a  source  of  fascination  for  phy- 
sicists ever  since  helium  gas  was  first  liquified  by  Kammer- 
lingh  Onnes  in  1908.   The  phase  diagram  of  helium  is  shown 
in  Fig.  1.   Above  the  transition  marked  the  X  line  in  Fig.  1, 
helium,  usually  referred  to  as  He  I,  despite  its  obvious 
quantum  properties,  shows  many  of  the  thermodynamic  and  hydro- 
dynamic  properties  characteristic  of  classical  liquids.   Be- 
low the  transition,  where  it  is  referred  to  as  He  II,  its 
characteristics  undergo  drastic  qualitative  changes  and  many 
of  its  properties  are  unique  to  helium.   The  most  outstanding 
of  these  changes  is  that  of  superfluidity — the  vanishing  of 
flow  resistance  for  the  passage  of  the  liquid  through  small 
channels.   These  properties  are  characteristic  of  natural 
helium,  which  is  nearly  100%   He.   The  lighter  isotope.   He, 
is  present  in  natural  helium  in  such  small  quantities  (about 
one  part  in  10  )  that  its  effects  are  not  noticeable  in  most 
experiments.   The  vast  amount  of  theoretical  and  experimen- 
tal work  on  liquid  helium  is  well  summarized  by  Wilks  (1967)  , 
by  Woods  and  Cowley  (1973) ,  and  by  Keller  (1969) . 

The  first  major  step  toward  a  theory  that  would  explain 
the  unique  properties  of  superfluid  helium  was  the  suggestion 
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Fig,  1.   Phase  diagram  of  ^He  [after  Wilks  (1967)]. 


made  by  London  (1954)  that  the  superfluid  phase  was  a  "fourth 
state  of  matter"  representing  a  macroscopic  manifestation  of 
quantum  effects.   He  reasoned  that  a  significant  fraction  of 

the  liquid,  or  perhaps  the  whole  liquid,  might  exist  in  a 

4 
single  quantum  state.   Because   He  has  zero  spin,  it  wxll 

obey  Bose-Einstein  statistics:  in  contrast,  the  rare  isotope 

3 

He  which  has  spin  1/2  will  obey  Fermi-Dirac  statistics.   In 

Bose  statistics  an  unlimited  number  of  particles  may  occupy 
an  energy  level.   This  has  the  mathematical  consequence  that 
an  ideal  gas  of  bosons,  i.e.,  a  large  number  of  non- 
interacting  particles  subject  to  Bose  statistics,  will  for 
certain  values  of  the  density  and  temperature  of  the  gas  have 
a  significant  fraction  of  the  particles  "condensed"  into  the 
ground  state.   This  "condensation"  will  occur  in  momentum 
space.   Thus  particles  separated  by  large  distances  in  the 
fluid  will  have  their  momenta,  hence  their  motion,  correlated. 

Fermi  statistics  do  not  permit  this  unless  there  is  a  pairing 

3 

of  the  atoms.   Thus  the  experimental  fact  that   He  is  not  a 

superfluid  until  less  than  0.003  K,  in  contrast  to  2.17  K  for 
the  superfluid  transition  in   He, is  easily  explained. 

This  suggestion  by  London  and  Tisza  provided  the  basis 
for  the  development  of  a  two-fluid  model.   London  (1954) 
gives  a  full  account  of  the  historic  development  of  this  idea 
and  the  role  played  by  Tisza  in  its  development.   Helium  was 
viewed  as  being  composed  of  two  interpenetrating  fluids--a 
superfluid  with  zero  entropy  and  a  normal  fluid.   The  two- 
fluid  model  successfully  predicted  most  of  the  macroscopic 


4 
thermal  phenomena  in   He.   It  also  explained  the  tremendous 

differences  in  low  temperature  properties  between   He  and 

4 
He.   However,  this  model  had  four  grave  weaknesses.   (1)  It 

was  a  macroscopic  and  not  a  microscopic  theory.   (2)  The 
agreement  between  theory  and  experiment  was  qualitative,  not 
quantitative.   (3)  The  theory  was  difficult  to  extend  to  dy- 
namical situations  such  as  the  propagation  of  sound.   (4) 
Finally,  the  theory  was  a  model,  many  of  whose  features,  such 
as  the  existence  of  the  condensate,  could  not  be  experimen- 
tally verified. 

The  next  major  advance  was  the  Landau  (1941,  19  47) 
theory  of  superfluidity.   Landau  viewed  superfluid  helium  as 
a  perfect  background  fluid  in  which  a  gas  of  elementary  exci- 
tations moves.   These  excitations,  which  have  a  definite 
energy  and  momentum,  describe  the  "collective  motion"  of  the 
helium  atoms  similarly  to  the  use  of  "normal  modes"  to  de- 
scribe the  motion  of  interacting  particles.   The  excitations 
behave  as  quasi-particles  which  interact  very  weakly  with 
each  other.   From  the  energy  spectrum,  i.e.,  the  relation  be- 
tween energy  and  momentum  for  the  excitations,  thermodynamic 
and  hydrodynamic  properties  of  liquid  helium  can  be  deduced. 
The  excitations  were  given  the  names  "phonons"  and  "rotons." 

Landau  pictured  the  phonon  excitations  as  long  wave- 
length density  fluctuations  with  energy,  e,    directly  propor- 
tional to  their  momentum,  p,  and  traveling  with  the  velocity 
of  sound,  c.   Thus, 


e  =  cp  .        (1.1) 

in  the  long  wavelength  limit.   At  shorter  wavelengths  (larger 
p)  the  e  versus  p  dispersion  curve  passes  through  a  minimum, 
such  that 

e  =  A  +  (p  -  p^)^/2y.  (1.2) 

The  excitations  in  the  region  of  the  minimum  were  called  ro- 
tons.   While  the  excitation  spectrum  is  a  smooth,  continuous 
curve,  at  low  temperatures  only  the  phonon  and  roton  states 
are  significantly  populated.   Since  the  excitations  were 
associated  with  the  "normal"  fluid  and  the  background  was 
identified  with  the  superfluid,  most  of  the  results  of  the 
two-fluid  model  could  be  derived  from  Landau's  theory.   How- 
ever, there  was  still  no  physical  microscopic  basis,  nor  was 
there  an  experimental  basis,  for  the  assumptions  of  the 
theory,  and  the  experimentally  verified  necessity  that  super- 
fluid  helium  be  a  Bose  system  had  vanished  from  the  theory. 
If  a  microscopic  physical  basis  for  the  energy  spectrum  of 
the  excitations  could  be  provided,  the  Landau  theory  would 
provide  an  adequate  description  for  almost  all  phenomena  in 
helium. 

Feynman  (1953,  1954)  shov/ed  that  a  liquid  consisting  of 
indistinguishable  particles  subject  to  Bose  statistics  could 
have  only  one  kind  of  low  energy  excitations.   These  were 
density  fluctuations  resembling  sound  v;aves  v;hich  v;hen  quan- 
tized were  identical  with  Landau's  phonons.   Since  Feynman 's 


arguments  were  very  general,  based  on  the  energy  of  a  con- 
figuration of  individual  atoms,  this  provided  a  microscopic 
explanation  for  the  phonon  part  of  the  excitation  spectrum. 
In  a  further  paper,  relying  again  on  very  generalized  argu- 
ments, Feynman  (1954b)  calculated  a  v/ave  function  for  liquid 
heD.ium  which  included  the  possibility  of  rotational  motion. 
By  applying  the  variational  principle,  the  energy  spectrum 
was  found  to  be 


2 

e(p)  =  |-  S(k)  (1.3) 


where  B(3:)  is  the  liquid  structure  factor.   For  the  small  k 
limit  at  T  -  0  this  reduces  to  e (p)  =  cp.   At  higher  momen- 
tum- the  structure  factor,  which  can  be  m.easured  by  X-ray 
scattering  from  the  fluid,  has  a  maximum  v/hich  leads  to  a 
minimum  in  e  (p) .   The  resulting  spectrum  is  in  qualitative 
agreement  with  Landau's  proposal.   An  improved  wave  function 
due  to  Feynman  and  Cohen  (1956)  greatly  improved  the  agree- 
ment between  the  computed  energy  spectrum  and  the  Landau  type- 
spectrum  as  determined  from  thermodynamic  measurements. 
Finally,  Feynman  and  Cohen  (1957)  pointed  out  that  the  exci- 
tations could  be  observed  directly  by  neutron  scattering. 
Neutron  scattering  experiments,  discussed  in  the  next  sec- 
tion, confirmed  the  excitation  picture  by  directly  measuring 
the  excitations  and  showed  that  their  spectrum  had  the  form 
suggested  by  Landau.   Thus  the  theory  had  a  firm  physical 


basis  and  the  condition  that  the  particles  obey  Bose  statis- 
tics was  automatically  included  in  the  theory. 

However,  this  was  still  not  a  "first  principles"  calcu- 
lation.  There  have  been  numerous  attempts  to  derive  the 
spectrum  from  a  consideration  of  the  interaction  potential 
for  two  helixim  atoms  and  the  application  of  statistical 
mechanics.   Despite  considerable  formal  progress,  the  micro- 
scopic theory  has  not  yet  been  developed  sufficiently  to  give 
a  detailed  description  of  liquid  he3ium.   However,  consider- 
able progress  has  been  made  in  variational  calculations  for 
the  structure  factor  and  the  ground  state  energy,  for  example, 
Pokrant  (1972)  .   A  comprehensive  siumtiary  of  these  efforts  and 
their  shortcomings  can  be  found  in  Keller  (1969)  and  in  Woods 
and  Cowley  (1972) . 

Neutron  and  X-ray  Scattex'ing  Measurements 
After  Feynman's  paper,  several  groups  performed  neutron 
scattering  experiments.   The  most  comprehensive  of  the  early 
works  was  that  of  Yarnell  et   at.     (1959) .   This  work  measured 
an  excitation  spectrum  that  was  in  good  agreement  with  the 
one  proposed  by  Landau  based  on  thermodynamic  arguments. 
Also  it  shov;ed  that  belov;  1.5  K  the  line  v;idth  of  the  exci- 
tations was  quite  narrow.   This  implied  that  the  excitations 
had  a  relatively  long  lifetime  and  that  the  interactions  be- 
tween excitations  were  weak,  two  critical  assumptions  of  the 
theory.   As  the  temperature  approaclied  the  X  line,  the  exci- 
tations became  very  broad  and  strongly  interacting  and 
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changed  character  for  temperatures  above  the  X  line,  thus 
confirming  that  the  measured  spectrum  was  connected  with 
superfluidity,  since  the  superfluid  properties  become  less 
pronounced  near  the  X    line  and  vanish  above  it. 

With  improvements  in  technique  and  technology,  more  work 
followed,  culminating  in  the  "definitive  work"  by  Cowley  and 
Woods  (CW,  1971)  ,  whose  v.-ork  was  accurate  enough  to  permit 
attempts  to  fit  the  spectrum  with  algebraic  expressions.   Re- 
cently Dietrich  et   at.     (1972)  undertook  a  comprehensive  study 
of  the  roton  dip  in  the  spectrum,  including  a  study  of  its 
variation  with  pressure.   Three  major  problems  have  occurred 
in  all  neutron  scattering  measurements.   Neutron  scattering 
is  limited  to  values  of  k  greater  than  about  0.4  A   ,  and  the 
measurem.ents  for  small  k  values  have  a  very  large  intrinsic 
error.   Except  for  the  roton  work  of  Dietrich  et   al.    and  some 
work  at  24  atm.  by  Svensson,  Woods  and  Martel  (1972) ,  and  Hen- 
shaw  and  Woods  (1961) ,  there  has  been  no  work  on  the  pressure 
variation  of  the  excitation  spectrum.   Finally,  all  neutron 
measurements  have  been  taken  above  1.1  K,  so  there  is  a  ques- 
tion of  possible  variation  of  the  curves  with  temper^iture  be- 
low that  point.   Nevertheless  neutron  scattering  is  the  most 
direct  and  most  accurate  way  of  measuring  the  excitation  spec- 
trum. 

In  addition  to  direct  measurement  of  the  excitation  spec- 
trum, neutron  scattering,  as  v/ell  as  X-ray  scattering,  can  be 
used  to  measure  the  structure  function,  S(k,a)).   The  structure 
factor  relates  the  differential  coherent  scattering  cross  sec- 
tion to  the  total  cross  section 


2  • 

o 


where  dQ  is  the  solid  angle  of  acceptance  of  the  scattered 

I 

beam,  dE  is  the  energy  width  of  the  scattered  beam, k  and  k 

are,  respectively,  the  incident  and  scattered  wave  numbers, 
o    is  the  total  cross  section  and  i^.o)  is  the  energy  transfer 
during  the  scattering.   Since  S  depends  on  the  momentum 
transfer  during  the  scattering  process,  frequently  S(k,a))  is 
written  S(Q,a))  where  Q  is  a  new  variable  defined  as  the 
momentum  transfer.   Because  the  Fourier  transform  of  S  is  the 
time  dependent  correlation  function,  S  is  related  to  the  spa- 
tial and  temporal  behavior  of  the  atoms  in  the  fluid. 

As  shown  by  Eq.  (1.3),  the  structure  factor  can  be  re- 
lated to  the  excitation  spectrum.   There  is  reasonable  agree- 
ment between  the  spectrum  obtained  directly  from  neutron 
scattering  and  that  deduced  indirectly  from  the  structure 
factor  as  obtained  in  the  best  recent  measurements  of  Hallock 
(1972) .   As  with  the  neutron  scattering  measurements  of  the 
excitation  spectrum,  structure  factor  measurements  are 
limited  to  k  greater  than  0.33  A~  .   There  has  been  essen- 
tially no  work  on  the  pressure  dependence  of  S.   Hallock  did 
carry  the  measurements  down  to  0.38  K,  which,  while  it  is  not 
T  =  0  as  used  in  most  theories,  does  eliminate  most  questions 
about  temperature  variation.   Scattering  can  also  be  used  to 
hunt  directly  for  the  presence  of  the  condensate,  but  results 
so  far  are  inconclusive  because  of  serious  experimental 
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problems--for  a  review  of  the  subject  see  Jackson  (1974) . 
Thus  scattering  measurements  have  thoroughly  demonstrated  the 
"reality"  of  Landau's  excitations  and  have  provided  a  quanti- 
tative basis  for  the  theory  which  is  in  good  numerical  agree- 
ment with  the  values  deduced  from  thermodynamic  measurements, 
but  there  are  still  gaps  in  the  measured  spectrum,  particu- 
larly in  the  long  wavelength-small  p  region. 

By  fitting  the  experimental  scattering  data  with  the 
assumed  form  given  by  Eqns.  (1.1)  and  (1.2),  explicit  values 
for  y.  A,  p   and  c  can  be  obtained.   These  values  can  then  be 
used  in  a  statistical  mechanical  analysis  to  yield  expres- 
sions for  the  pressure,  specific  heat,  entropy,  expansion  co- 
efficient and  the  superfluid  density  used  in  the  two-fluid 
model.   This  will  be  done  explicitly  in  Chapter  II.   These 
expressions  can  be  compared  directly  with  experimental 
values.   Historically,  the  process  was  reversed,  with  values 
for  y,  A,  p   and  c  being  deduced  from  thermodynamic  measure- 

0 

ments  and  these  values  were  used  to  construct  the  excitation 
spectrum.   In  general  there  is  good  agreement  between  calcu- 
lated and  measured  thermodynamic  properties.   This  agreement 
is  usually  expressed  by  comparing  the  values  of  y,  p  ,  A  and 
c  derived  from  scattering  with  those  derived  from  thermo- 
dynamics.  There  is  a  major  flaw  in  this  scheme:   since 
thermodynamic  measurements  simultaneously  sample  the  entire 
excitation  spectrum,  any  deviation  of  the  spectrum  from  the 
assumed  forms  can  lead  to  a  change  in  the  parameters.   This 
is  a  particularly  severe  problem  with  rotons,  since  the 
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measured  excitation  spectrum  is  not  parabolic  except  in  a 
very  narrow  region  around  the  minimum.   Most  of  the  differ- 
ence between  the  thermodynamic  values  for  \i ,    A,  p   and  c,  and 

0 

the  neutron  values  of  y,  A,  p   and  c  can  be  accounted  for  by 

0 

this  problem.   Thus  the  phonon-roton  picture  of  liquid  helium 
rests  on  a  sound  phenomenological  basis — neutron  scattering — 
and  successfully  predicts  most  of  the  thermodynamic  proper- 
ties of  helium. 

Measurements  of  the  Velocity  and 
Attenuation  of  Sound 

Only  one  parameter  in  the  phonon-roton  model  can  be 

directly  measured  other  than  by  scattering.   That  parameter 

is  c,  the  sound  velocity.   Many  measurements  of  the  sound 

velocity  and  the  attenuation  of  sound  have  been  made.   The 

most  comprehensive  of  these  was  a  series  of  measurements  made 

at  the  Argonne  Laboratories,  Abraham  et   at.     (1969,  1970)  and 

Roacli  ct   al.     (1972a,  1972b)  ,  which  measured  sound  velocity, 

attenuation  of  sound,  and  the  Gruniessen  constant,  V   =   —   ^r— 

C  dp 

(where  p  is  the  density) ,  as  a  function  of  pressure,  fre- 
quency and  temperature  in  the  region  below  0.6  K.   These 
measurements  could  be  compared  V7ith  the  detailed  predictions, 
for  these  quantities  in  liquid  helium,  developed  by  Khalat- 
hikov  and  co-workers — Andreev  and  Khalatnikov  (1963) ,  Khalat- 
nikov  (1965),  and  Khalatnikov  and  Chernikova  (1965,  1966). 
Their  results  were  based  on  the  development  of  quantum  hydro- 
dynamics for  superfluid  helium  by  Landau  and  Khalatnikov  as 
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summarized  by  Khalatnikov  (19  65) .   For  there  to  be  attenu- 
ation of  sound,  Eq.  (1.1)  was  modified  to 

e  =  cp(l  -  Yp^).  (1.5) 

Some  modification  was  necessary,  since  a  pure  phonon  spec- 
trum, e  =  cp,  does  not  permit  phonon  decay  because  of  the 
impossibility  of  simultaneously  satisfying  energy  and  momen- 
tum conservation  except  through  the  quantum  mechanical  un- 
certainty of  c.   Experimentally  the  energy  width  is  known  to 
be  far  too  small  to  permit  the  observed  attenuation.   The 
additional  term  containing  y  is  known  as  dispersion.   The 
exact  form  of  the  dispersion  term  was  a  reasonable  assump- 
tion, but  was  not  the  only  possibility.   Since  Landau  and 
Khalatnikov' s  theory  had  been  successful  in  predicting  most 
hydrodynamic  phenomena  in  liquid  helium,  it  was  a  major  sur- 
prise when  their  predictions  were  qualitatively  wrong  about 
the  measured  values  of  the  velocity  and  attenuation  of  sound 
as  well  as  their  variation  with  temperature,  pressure,  and 
frequency  at  very  low  temperatures.   The  theoretical  predic- 
tion for  the  attenuation,  a,  was  that  in  low  temperature 

limit  a  would  be  proportional  to  T  .   The  experimental  result 

4 
was  that  a  was  proportional  to  T  .   Furthermore  the  measured 

attenuation  was  approximately  two  times  greater  than  the 

Calculated  attenuation.   As  the  pressure  was  raised,  a 

strange  "shoulder"  appeared  in  the  data  although  nothing  of 

the  sort  had  been  predicted.   The  velocity  of  sound  showed 

extremely  complicated  behavior  as  a  function  of  frequency  and 
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temperature  which  had  not  been  predicted.  As  Abraham  et  al. 
(1969)  remarked,  "We  therefore  conclude  that  the  present 
theoretical  formulation  of  sound  propagation  at  very  low 
temperatures  is  incomplete"  (p.  370)  .  How  could  a  model 
which  had  so  successfully  predicted  so  many  of  the  properties 
of  helium  fail  so  miserably  in  the  lov;  temperature  limit 
where  it  should  perform  the  best? 

Recent  Theoretical  Developments 
Although  substantially  confirmed  by  neutron  scattering, 
the  phonon-roton  picture  of  helium,  like  the  two-fluid  model, 
is  a  model.   Keller's  (1969)  warning  about  models  is  appro- 
priate.  "A  model  is  devised  to  represent  a  physical  phenom.e- 
non  because  the  actual  situation  is  far  too  complicated  to  be 
handled  directly;  the  model  then  incorporates  the  simplifying 
assumptions  that  make  the  problem  tractable.   The  first  dan- 
ger is  one  of  oversimplification,  and  this  is  an  especially 
hazardous  possibility  when  quantum  effects  are  involved  .... 
When  such  a  [good]  model  exists,  there  is  a  tendency  to 
take  it  too  literally,  to  promote  it  to  too  high  a  status, 
and  then  finally  to  forget  it  is  a  model"  (p.  17) .   One  such 
oversimplification  of  the  phonon-roton  model  can  be  seen  in 
Fig.  2.   Clearly  phonons  and  rotons  are  a  very  poor  approxi- 
mation to  the  actual  excitation  spectrum  in  the  region  of  the 
maximum  of  the  curve  and  on  the  high  k  side  of  the  minimum. 
To  check  the  possibility  that  this  oversimplification  might 
affect  calculations  based  on  the  model,  Bendt  et   al.     (1959) 
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Fig.  2.   Excitation  spectrum  at  SVP  [after 
Cowley  and  Woods] . 

dashed  line  is  the  phonon  approximation, 
dashed  and  dotted  line  is  the  roton 
approximation . 
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and  Singh  (1968)  made  the  laborious  and  complex  calculations 
of  some  of  the  thermodynamic  properties  of  liquid  helium 
directly  from  the  measured  excitation  spectrum.   They  found 
that  calculations  based  on  the  phonon-roton  model  deviate 
very  little  from  their  exact  calculations  at  low  tempera- 
tures, but  that  those  deviations  increase  sharply  with  in- 
creasing temperature,  becoming  quite  large  above  1.5  K.   Be- 
cause the  v;hole  model  collapses  rapidly  as  one  approaches  the 
X  transition,  this  is  not  surprising.   Unfortunately,  no  work 
has  ever  been  donc>  on  the  pressure  dependence  of  the  devia- 
tions of  the  phonon-roton  model  from  exact  calculations, 
partly  because  no  extensive  neutron  scattering  studies  have 
been  done  as  a  function  of  pressure.   Only  very  high  preci- 
sion measurements  of  thermodynamic  properties  will  show  de- 
viations from  the  phonon-roton  model  since  the  total  devia- 
tion is  of  the  order  of  2%  for  temperatures  near  1  K  and  de- 
creases rapidly  as  the  temperature  is  lowered  below  1  K. 

There  would  appear  to  be  little  possibility  that  there 
is  any  major  flaw  in  the  theory  of  sound  propagation,  since 
it  has  v/orked  so  v;ell  for  other  materials.   Thus  attention 
focused  on  the  approxim.ation  that  e  (p)  =  cp(l  -  YP")  •   Iff  as 
Landau  had  assum.ed,  y   v;as  greater  than  0,  then  the  dispersion 
was  normal,  i.e.,  the  type  found  in  most  materials.   Normal 
dispersion  prohibits  three-phonon  processes  by  the  conserva- 
tion laws.   Three-phonon  processes  are  those  where  one  phonon 
decays  into  two  others,  or  v;here  tv70  phonons  combine  to  make 
a  third.   All  theorists  agree  that  the  only  way  to  produce 
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4 
the  experimentally  observed  T  dependence  of  the  attenuation 

of  sound  at  low  temperatures  is  by  allowing  three-phonon  pro- 
cesses to  occur. 

Simon  (1963)  and  Petick  and  Ter  Haar  (1966)  showed  that, 
under  some  circumstances,  three-phonon  processes  were  allowed 
without  dispersion.   Due  to  the  finite  lifetime  of  the  pho- 
nons,  the  phonon  spectrum  has  a  width  which  allows  the  con- 
servation equations  to  be  satisfied.   However,  Friedlander, 
Eckstein  and  Kuyper  (1972)  claimed  to  have  shown  that  proper 
renorraalization  of  phonons  forbids  three-phonon  processes  for 
all  cases  of  "normal"  dispersion.   The  best  fit  to  the  CW  re- 
rults  at  saturated  vapor  pressure  was  the  phonon  spectrum 
e  =  cp(l  -  Yp^  -  &p^)    with  Y  =  0  ±  2  X  10^^  and 
5  =  2.4  ±  .3  X  lo"^^.   This  disagreed  sharply  with  y  =  8  x  10^^ 
derived  from  ultrasonic  measurements  by  Eckstein  and  Varga 
(1968)  using  Khalatnikov' s  approach.   Havlin  and  Luban  (1972) 
were  able  to  fit  the  ultrasonic  data  with  a  spectrum 
e (p)  =  cp(l  -  6p^);  but  their  6  -  0.9  x  10^^  was  a  factor  of  3 
lower  than  the  neutron  value.   As  will  be  shown  in  detail  in 

Chapter  II,  no  thermodynamic  measurement  can  separate  the 

2 

effect  of  Y   from  the  effect  of  5.   So,  unless  y  remains  very 

small  over  a  considerable  range  of  pressure,  an  analysis  based 
solely  on  5  has  limited  usefulness.   The  high  pressure  neu- 
tron measurements  of  Svensson,  Woods  and  Martel,  which  had 

37 
Y=6.2±0.6xl0   ,  argued  against  any  assumption  that  y 

was  insignificantly  small  over  the  entire  pressure  range  of 

liquid  helium. 
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Maris  and  Massey  (1970)  suggested  that  if  y  were  less 
than  zero — so-called  anomalous  dispersion--the  problems  of 
the  theory  would  be  solved.   Anomalous  dispersion  unques- 
tionably permits  three-phonon  processes.   Further  calcula- 
tions by  Maris  (1372,  1973),  in  which  he  solved  the  integral 
equations  for  velocity  and  attenuation  of  sound  numerically 
without  the  customary  approximations,  but  assumed  anomalous 
dispersion,  were  in  good  agreement  with  the  experimental  data 

on  attenuation  and  velocity  of  sound  below  about  14  atm. 

37 
pressure.   These  calculations  yielded  a  value  of  -8  x  10 

for  Y  ^t  saturated  vapor  pressure. 

Jackie  and  Kehr  (1971)  v/ere  able  to  explain  the  shoulder 
in  the  attenuation  data  by  using  anomalous  dispersion  and  in- 
troducing a  cutoff  frequency.   The  cutoff  occurred  because  of 
an  assumption  that  an  ultrasonic  phonon,  i.e.,  one  produced 
by  the  experimenter,  could  only  be  absorbed  by  a  thermal  pho- 
non, i.e.,  one  already  present  in  the  liquid  due  to  its  fi- 
nite temperature,  if  the  thermal  phonon  had  a  momentum  "Rk 
less  than  "hk^,.   This  reproduced  the  shoulder  very  well  if  k^ 
was  a  strong  function  of  pressure,  decreasing  sharply  with 
increasing  pressure.   That  implied  that  the  shape  of  the  pho- 
non spectrum  was  very  sensitive  to  pressure.   Using  available 
data  on  the  shape  of  the  curve,  they  predicted  that  the 
effect  of  the  cutoff  would  be  most  noticeable  between  14  and 
19  atm. ,  in  accord  with  experiment. 

The  assumption  that  y  is  negative  would  radically  change 
the  generally  accepted  ideas  about  phonon-phonon  interactions. 
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V^Tiile  agreement  with  the  ultrasonic  data  was  gratifying,  at 
that  time  there  was  little  other  experimental  evidence  to 
support  this  assumption.   There  were  hints  from  the  neutron 
data  that  this  might  be  possible,  since  earlier  neutron  work 
by  Henshaw  and  VJoods  (1961)  had  reported  the  possibility  that 
Y  v/as  negative.   There  was  also  evidence  from  X-ray  scatter- 
ing that  the  small  k  structure  factor  was  not  inconsistent 
with  a  negative  y.   This  result  was  given  a  boost  with  a 
model  calculation  by  lachello  and  Rassetti  (1973)  ,  based  on  a 
new  technique  of  deriving  the  spectra  from  the  helium  poten- 
tial curve,  yielded  a  model  which  had  anomalous  dispersion  at 
low  pressures  and  normal  dispersion  at  higher  pressures. 
Also  Zasada  and  Pathria  (1972)  had  shown  a  similar  result; 
negative  y  at  low  densities,  positive  y  at  higher  densities 
could  be  obtained  for  an  imperfect  Bose  gas  by  using  reason- 
able potentials.   Experimental  evidence  discussed  in  the  next 
subsection  has  tended  to  support  this  hypothesis. 

A  recent  calculation  by  Jackie  and  Kehr  (1974)  concluded, 
on  the  basis  of  ultrasonic  measurements  at  finite  tempera- 
tures, that  y  was  negative  and  nearly  constant  v/ith  pressure, 
while  6  was  positive  and  increased  rapidly  with  pressure. 
Thus  the  apparent  normal  dispersion  at  higher  pressures  would 
be  caused  by  5  becoming  so  large  as  to  dominate  y  rather  than 
by  y  changing  signs. 

All  of  the  preceding  calculations  have  assumed  that  the 
phonon  spectrum  is  given  by 
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e(p)  =  cpd  -  Yp^  -  5p^...).  (1.6) 


The  series  is  ass\imed  to  terminate  with  the  5  term  in  the 
sense  that  higher  order  terms  are  too  small  to  be  observed. 
Feenberg  (1971)  showed  that  for  a  simple  interatomic  poten- 
tial, such  as  the  inverse  r  ,  e.g.,  the  dipole-dipole  poten- 
tial, that  S (p)  and  e (p)  must  contain  both  even  and  odd 
pov/ers  of  p  beyond  the  cubic  term.   Thus  there  was  no  a 
priori  reason  why  the  expansion  for  e (p)  should  not  include 
all  pov/ers  of  p.   This  suggestion  led  Molinari  and  Regge  (MR, 
1971)  to  reanalyze  the  results  of  CVl   yielding  a  five-term  fit 


e(p)  =  cmpd  +  0.5465p  -  1.3529p^  +  0.2595p^  +  O.lSSOp'^ 

-  0.0522p^)^/^  (1.7) 

with  both  even  and  odd  powers.   Most  experimental  evidence 

2 

seems  to  indicate  that  the  coefficient  of  the  p   term  of  the 

expansion  is  zero.   This  destroys  the  ability  of  the  formula 
to  fit  the  neutron  data. 

At  about  the  same  time,  Gould  and  Wong  (1971)  showed 
that  for  a  weakly  interacting  Bose  gas  there  are  terms  of  the 
form  p  log  (1/p)  in  both  e (p)  and  S  (p) .   Their  fairly  general 
arguments  indicated  that,  to  the  extent  that  helium  could  be 
treated  as  a  weakly  interacting  Bose  gas,  no  analytic  expan- 
sion of  e (p)  was  possible.   However,  since  their  complete 
spectrum  was  of  the  form  e  (p)  =  CqP  +  C2P   +  c^^p   log  (1/p)  it 
would  be  difficult  to  detect  the  non-polynomial  terms.   The 
source  of  the  singularity  lay  in  the  multiphonon  interactions. 
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While  multiphonon  interactions  are  important  in  helium,  it  is 
very  questionable  whether  a  weakly  interacting  Bose  gas  is  a 
good  model  for  multiphonon  interactions  in  helium.   So  far, 
no  one  has  analyzed  the  experimental  data  to  look  for  evidence 
of  the  predicted  log  (1/p)  term,  due  to  the  formidable  mathe- 
matical difficulties  in  doing  such  an  analysis. 

Lin-Liu  and  Woo  (1974)  presented  a  calculation  based  on 
sum  rules  which  the  structure  should  obey.   This  yielded  a 
phonon  spectrum  at  P  =  0  and  T  =  0  of  the  form 

e(p)  =  cp(l  +  0.17  ^^   +   0.78  -^  ^  -  3.3  ^-  ^).  (1.8) 

^     '^  n\c  mc         mc 

The  positive  numbers  indicate  anomalous  dispersion.   In  this 

3 

calculation,  the  anomalous  dispersion  comes  from  both  the  p 

term — the  term  with  coefficient  y  in  the  usual  expansion — and 

4  5 

the  p   term.   The  p   term — the  term  with  coefficient  6  m  the 

usual  expansion — has  normal  dispersion.   No  convments  were 

made  about  the  pressure  dependence  of  any  of  the  coefficients. 

With  the  exception  of  the  work  of  Jackie  and  Kehr  (1974) ,  all 

calculations  of  the  spectra  are  for  T  =  0  at  zero  applied 

pressure.   With  the  exception  of  MR,  all  of  the  formulas  for 

the  excitation  spectrum  are  in  agreement  v/ith  experimental 

data.   This  is  hardly  surprising  since  most  experiments  can 

2 

just  barely  detect  the  effect  of  the  yp   term,  making  posi- 
tive identification  of  the  higher  order  terms  difficult  and 
the  assignment  of  accurate  nuinerical  values  virtually  impos- 
sible. 
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Recent  Experimental  Results 
The  principal  experimental  evidence  has  been  a  set  of 
specific  heat  measurements  by  Phillips,  Waterfield  and  Hoffer 
(PWH,  1970) .   They  measured  the  specific  heat  of  liquid 

helium  at  four  different  pressures.   At  low  pressures  they 

2 
found  evidence  for  anomalous  dispersion.   If  e  =  cp (1  -  yp  ) , 

37 
then  Y  =  -  4.1  X  10    at  saturated  vapor  pressure.   y  de- 
creased in  magnitude  with  increasing  pressure  and  became 

positive  somewhere  between  5  and  20  atm.   At  20  atm.,  they 

•37 
found  Y  =  +  19-6  X  10   .   They  claimed  their  work  could  also 

be  fitted,  though  "less  well"  at  saturated  vapor  pressure  by 

7  fi 

Y  =  0  and  6=-4.5xl0.  No  other  fits  were  attempted. 
There  is  some  question  about  the  validity  of  the  analysis 
performed  on  the  specific  heat  data  by  PWH.  Nevertheless, 
the  data  are  so  significant  that  most  theories  are  tested  by 
reanalyzing  these  data  to  see  if  they  can  be  fitted  by  the 
theory.  No  other  thermodynamic  data  accurate  enough  to  be 
analyzed  have  been  published. 

A  recent  reanalysis  of  these  specific  heat  data  by 
Zasada  and  Pathria  (1974) ,  using  the  measured  neutron  spec- 
tra instead  of  the  phonon-roton  approximation,  yielded 
37 

Y  =  -  5.1  X  10    at  saturated  vapor  pressure.   Similar 
analysis  was  performed  on  the  data  at  the  other  two  "low 
pressures,"  changing  y  slightly,  but  no  fit  v/as  attempted 
for  the  highest  pressure. 

Similarly  the  X-ray  measurements  of  S (p)  of  Hallock 

37 
(1972)  could  be  analyzed  to  give  y   =   -    5.1   x   10      .      However, 
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they  could  also  be  analyzed  to  yield  a  quadratic  term  as  pro- 
posed by  MR  or  to  give  aY=-  3x10   ,  i.e.,  =0  to  within 

the  quoted  error  of  CW. 

2 
Since  the  quadratic  term  a,p  ,  proposed  by  MR,  was  of 

3 

lower  order  than  the  usual  yp   term,  there  were  immediate 

attempts  to  observe  it.   Anderson  and  Sabisky  (1972)  measured 
the  acoustic  thickness  of  helium  films  at  1.38  K  for  frequen- 
cies between  20  and  60  GHz.   Their  analysis  yielded  a  posi- 
tive quadratic  term  in  good  agreement  with  MR.   These  frequen- 
cies were  over  an  order  of  magniture  greater  than  the  usual 
ultrasonic  frequencies,  the  largest  of  which  was  256  MHz. 
The  temperature  was  also  much  higher  than  usual,  since  most 
measurement  of  phonon  properties  are  made  below  0.6  K.   Re- 
cently Anderson  and  Sabisky  (1974)  stated  that  their  analysis 
of  film  properties  may  have  been  inadequate. 

Roach  et   al.     (1972c)  measured  the  frequence  dependence  of 
the  sound  velocity  with  an  improved  technique  at  T  =  0.3  K  at 
frequency  30  and  9  0  MHz.   They  reported  that  if 
£  (p)  =  cp(l  +  a,p) ,  then  a,  =  0  ±  0.01  A,  in  sharp  contrast  to 
a^  =  0.275  ±  0.030  A  for  Anderson  and  Sabisky  (1972)  and  MR. 
A  reanalysis  of  the  specific  heat  data  by  Zasada  and  Pathria 
(1972)  concluded  that  the  inclusion  of  a  quadratic  term  in 
the  energy  spectrum  led  to  a  set  of  parameters  whose  behavior 
as  a  function  of  density  was  too  erratic  to  be  acceptable; 
whereas,  setting  a,  identically  equal  to  zero  led  to  a  set  of 
parameters  which,  as  a  function  of  density,  v;ere  in  agreement 
with  direct  measurements  of  the  same  quantities.   This 
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experimental  evidence  is  generally  regarded  as  showing  that 
there  is  no  quadratic  term  in  the  excitation  spectrum. 

The  uncertainties  in  the  experimental  situation  are 
illustrated  by  two  recent  experiments.   (1)   Narayanamurti, 
Andres  and  Dynes  (1973)  measured  the  group  velocity  and 
attenuation  of  phonons  with  frequencies  of  2  x  10    to  9  x 
10"   Hz,  far  in  excess  of  any  ultrasonic  measurements.   They 
report  zero  dispersion  over  the  whole  range  of  frequencies. 
As  noted  by  Jackie  and  Kehr  (1974) ,  there  appear  to  be  some 
inconsistencies  in  the  analysis  of  Narayanamurti  et   al.       For 
example,  in  1974  Dynes  and  Narayanamurti  used  three-phonon 
processes  to  analyze  their  experiments  with  heat  pulses  in 
helium.   The  phonon  lifetimes  used  in  the  heat  pulse  analysis 
are  inconsistent  v.'ith  three-phonon  processes  occurring  in  the 
absence  of  dispersion.   The  heat  pulse  experiments  did  show  a 
striking  change  near  15  to  17  atm.  and  normal  dispersion  at 
higher  pressures.   (2)   Mills  et   al.     (1974)  measured  the 
angular  spreading  of  phonon  beams,  which  they  interpreted  to 
be  evidence  for  negative  (anomalous)  dispersion  balow  17  atm. 
and  normal  dispersion  above  that.   This  analysis  was  based  on 
the  presence  or  absence  of  three-phonon  processes  as  shown  by 
the  width  of  a  phonon  beam. 

In  summary,  the  weight  of  the  evidence  suggests  that 
there  is  anomalous  dispersion  in  liquid  helium  at  low  pres- 
sures and  normal  dispersion  at  higher  pressures.   However, 
more  high  precision  experimental  data  are  needed  to  decide 
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among  the  different  theories.   This  thesis  was  undertaken  to 
provide  high  precision  measurements  of  the  thermodynamic 
pressure  as  a  function  of  both  density  and  temperature  below 
1  K.   The  results  will  be  compared  with  the  predictions  of 
the  different  theories.   The  measurements  were  also  carried 
to  higher  temperatures  to  compare  them  with  the  measurements 
of  pressure  as  a  function  of  density  and  temperature  above 
1.5  K  made  by  Keesom  and  Keesom  in  19  33. 


CHAPTER  II 
LANDAU  THEORY 

Algebraic  Results 
In  this  chapter  the  predictions  of  the  thermodynamic 
pressure  for  the  different  phonon  spectra  will  be  computed. 
If  one  adopts  the  Landau  model  for  helium,  then  the  thermo- 
dynamic properties  of  liquid  helium  are  determined  by  the 
excitations.   If  there  are  Np  excitations  with  m.omentum  p 
and  energy  e  (p) ,  the  total  energy  of  the  fluid  is 


E  {N^}  =  E^  +  E  N^  e(p)  (2.1) 

P      o    p   p    ~ 


where  E   is  the  ground  state  energy.   This  is  the  fundamental 
assumption  of  Landau--all  of  the  energy  of  the  system  can  be 
accounted  for  by  the  excitations  of  the  single  particle  exci- 
tation spectrum  e  (p) .   However,  this  is  not  strictly  true. 
As  CW  showed,  there  are  other  branches  of  the  excitation 
spectrum.   However,  these  lie  at  such  high  energy,  e/k  greater 
than  18  K,  that  they  should  have  little  effect  on  low  tempera- 
ture (less  than  1  K)  properties,  due  to  the  v/eighting  factor 
exp  (-e/kT)  in  thermodynamic  calculations. 

By  definition,  the  partition  function  for  a  liquid  is 


Z  =  Z 


„  E   exp(-3E{N„})  (2.2) 

{Np}  ^ 
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where  3  has  its  usual  definition,  3  =  1/kT.   Assuming  that 
the  excitations  are  bosons,  v;hich  would  be  expected  for  a 
system  of  Bose  particles,  this  equation  reduces  to 

Z  =  Z__n[l-exp(-Be{p))]"^  (2.3) 

Op 

The  Helmholtz  free  energy  is  defined  to  be  F  =  -  kT  log  Z. 
Thus,  F  is  given  by 

F  =  -kT  Z  log[l-exp(-3e  (p)  )  ]~"^  +  Fq  (2.4) 

P 

where  Fq = -kT  log  Z^  is  the  ground  state  free  energy. 

Assuming  a  continuum  distribution  of  states  so  that  T.   - 

P 

-^  /  d^p  ,  then 


F  =  -  ~f  /  log[l-exp(-Be(p))]  -"d-^p  +  F^         (2.5) 


This  result  assumes,  in  agreement  with  experimental  results, 
that  e  (p)  is  a  smooth,  v;ell-behaved  function.   For  a  fluid 
one  can  assume  isotropy,  i.e.,  properties  of  the  material  de- 
pend only  on  the  magnitude  of  p,  not  its  direction.   Thus  one 

3        2 
can  replace  d  p  by  47Tp  dp,  yielding 


47rkTV   °°  7 

F  =  Fo  +    \        r  log[l-exp(-3c(p))]p^dp.  (2.6) 


h 


3 


Integrating  by  part?  gives 


3 

F  =  F   -  ^^'^      5      P   d[Bc(p)]  .„  ... 

o     3h3    oexp(3e"(p))-l-  ^^.Z) 
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This  assumes  that  Lim  ^e (p)  =  ^,    which  is  true  for  all  phy- 
sical systems. 

In  addition  to  the  mathematical  approximations  discussed 
above,  there  are  two  assumptions  inherent  in  this  approach 
which  must  have  their  validity  checked.   (1)   The  excitations 
are  well-defined,  long-lived  and  weakly  interacting  so  that 
they  can  be  treated  as  independent  entities.   This  is  con- 
firmed by  neutron  measurements.   (2)   The  excitation  spectrum 
is  a  fixed  quantity,  neither  varying  with  temperature  for  a 
fixed  density  nor  varying  rapidly  with  small  density  changes. 
To  be  useful,  the  spectrum  must  not  vary  during  an  experi- 
ment.  For  example,  it  v;ill  be  assumed  that  the  spectrum  is 
constant  in  this  experiment,  although  the  temperature  varies 
from  0.3  K  to  0.9  K.   Since  neutron  scattering  measurements 
extend  only  to  1.1  K,  this  assumption  cannot  be  experimen- 
tally proven  for  temperatures  less  than  1  K.   However,  neu- 
tron measurements  are  almost  independent  of  temperature  be- 
tween 1  K  and  1.5  K.   Similarly,  the  sound  velocity,  which  is 
an  explicit  parameter  in  the  excitation  spectrum,  varies  by 
less  than  0.1%  for  temperatures  between  0.1  K  and  1  K  and 
varies  slowly  with  density.   Hovrever,  theoretical  calcula- 
tions by  Ishikawa  and  Yamada  (197  2)  showed  that  for  an  im- 

2 
perfect  Bose  gas  the  higher  order  terms,  such  as  yp  /  varied 

considerably  with  temperature  and  that  temperature  variations 

in  the  interactions  between  excitations  v,'ould  have  the  saiae 

effect  as  introducing  new  higher  order  terms.   Because  these 

variations  only  become  large  for  temperatures  greater  than  1 
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K,  the  rest  of  this  development  of  the  theory  will  assume 
that  the  excitation  spectrum  is  fixed  in  temperature. 

Knowing  e (p)  from  neutron  and  X-ray  scattering,  and  more 
indirect  sources,  one  calculates  F  from  Eqns •  (2.6)  or  (2.7). 
Then,  the  thermodynamic  functions  P,  S,  E,  and  C^  are  calcu- 
lated by  the  usual  formulas: 

P  =  -   (|f,^,   S  =  -   (|f,^,   E  .  F  +  TS,  C^  =   (||)^. 

In  practice,  the  calculation  of  F  is  very  difficult  unless 
further  approximations  are  made.   Furthermore,  the  calcula- 
tion of  the  thermodynamic  functions  is  extremely  difficult 
unless  F  is  known  as  an  algebraic  function  of  temperature, 
T,  and  volume,  V. 

The  usual  approximation  for  helium  at  saturated  vapor 
pressure  is  shown  by  the  dotted  lines  in  Fig.  2.   Fig.  3 
shov/s  the  actual  excitation  spectrum  and  the  phonon-roton 
approximation  for  helium  at  2  4  atm.  pressure.   Separate  free 
energies  are  computed  for  the  phonons  and  rotons.   Thus,  F  = 
^o  ■*■  ^r  +  ^ph*   ^O'  the  free  energy  at  0  K,  is  independent  of 
temperature.   The  volume  derivative  of  Fq  gives  Pq,  the  pres- 
sure at  0  K.   Pq  will  be  treated  as  an  experimental  parameter 

and  will  not  be  calculated.   To  calculate  Fj-  one  inserts  the 

2 
roton  spectrum  e  (p)  =  A  +  (p-Po)  /2y  into  Eq.  (2.6).   Only  a 

very  small  error  is  introduced  by  extending  the  limits  of 

integration  to  0  and  0°.   Since  the  energy  gap,  A,  is  large 

compared  to  the  temperature,  the  logarithm  can  be  expanded  to 
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first  order,  i.e.,  log  [1-exp  (-ge (p) )  =  exp (-Be (p) ) ] .   Because 
of  the  exponential  weighting  factor  and  the  relatively  small 

width  of  the  minimum,  the  major  contribution  to  the  integral 

2 

will  come  from  p'\'p  ,  so  that  the  p   m  the  integrand  of  Eq. 

2 
(2.6)  may  be  replaced  by  p   .   This  gives 

F^  =  -  ^1}^   exp(-AAT)p^^  /"  exp(-(p-p^)V2ykT)dp.      (2.8) 


Integrating  gives  the  result 

Fj.  =  -(2(ykT)-'-/^p^^kTV  exp  (-A/kT)  )  /  (  (2ir)  ^/"Ti"^)  .       (2.9) 

This  yields 
-9F 


P   = 


8V 


-  =  2  (ykT)-'-'^^p^kTV  exp(-A/kT)  (2.10) 


U  +  A  £  |A       P  |M  _  2p  ^Po3/((  ^^3/2^.3) 
kTA3p   ^   y3p   p   9p 


In  order  to  compute  the  pressure  in  a  useful  form,  the 
variation  of  A,  p   and  y  with  density  must  be  known.   Donnelly 

(1972)  calculated  the  following  values  for  the  Landau  parame- 

3 

ters  as  a  function  of  density.   For  p  in  gm/cm  ,  these  values 

are 


A/k  -  (16.99  --  57.31p)  Kelvin 

1/3 
p  /fi  =  3.64p  ^   inverse  Angstroms  (2.11) 

y  =  (0.32  -  1.103p)M/    grains 

^He 

These  simple  algebraic  forms,  based  on  the  neutron  scattering 
work  of  Dietrich  et   al.     (1972)  will  be  used  in  preference  to 
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the  graphical  tables  of  Mills  (1965)  ,  which  v/ere  baaed  on 
thermodynamic  data.   At  loV'/  pressure  the  two  different 
raethods  yield  similar  values  for  A,  p^  and  ;i ,  but  have  sub- 

sbantially  different  values  for  the  density  derivatives  ^ — / 
3p        jv 

-j5 — ,  and  -— .   At  higher  pressures,  the  two  approaches  dis- 
agree about  both  the  values  of  the  Landau  parameters  and  the 
values  of  their  density  derivatives.   The  differences  between 
the  Mills  and  the  Donnelly  roton  spectra  can  be  seen  in  Fig. 
3.   Measurements  of  Sj-  by  van  den  Meijdenberg  eb   at.     (1961) 
yield  Landau  parameters  in  substantial  agreement  with 
Donnelly  if  the  corrections  at  low  pressures  suggested  by  van 
den  Meijdenberg  et   at.    are  applied  to  their  results.   Using 
Donnelly's  values  for  A,  Pq  and  y  and  their  variation  with 
density,  one  obtains 

P   =    ^^'^^  ..,  {[kT(0.32  -  1.103p)]^/^  [3.64p^/^]^       (2.12) 
^     (2iT)3/2n 

[1/3  -  ^^-^  +  ^_^4515p__^  fexp((57.31p  -  16.99)/T)]} 

Fpj^  is  calculated  in  a  similar  fashion.   A  phonon  spectrum, 
the  most  general  form  encompassing  all  proposed  formulas  ex- 
cept that  of  Gould  and  Wong, 

2       3       4 
e  (p)  =  cp(l  +  a,p  +  a2P   +  a3p   +  a^p  ...) 

is  inserted  in  Eq.  (2.7).   The  limits  of  integration  are 
extended  to  infinity.   This  causes  complete  overlap  of  the 
integrals  for  Fj.  and  F^.   However,  the  resulting  error  is 
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small  because  of  the  large  values  of  e   and  e  ,  in  the  re- 

r      pn 

gions  of  overlap.   This  yields 

ph  To/-iv,33  -5  c  7  12c 

*^  180h   c  27T 

^    15B(7),,c  1  Tc       3^/2uk>^_7 

+   ---y-(15a^a2    -    3a3    -    16a^    )  (--)    T 

Stt  ' 
+    l/2(54aj^^    +    12a2^    -    75aj^^    +    ISa^a^    -    3a^) 

(^^)    T*^]  (2.13) 

where  B(n)  equals  the  Gamma  function  of  n  times  Zeta  function 
of  n.   Hence,  the  pressure  is 

Pnh=--T4^(l+3r)T^    -   i^^    (2lll)[(l  +  4r)a,  -p   ^IT^ 
P^      ISOh^c^  2ti5  ^  1  3P 

+    10/7     (^)^     [(1  +  51)  {3a3^^    -a2)  -  P(6a3^   -y^--y^)]T^ 

+   -^r  B(7)  (^^)        [(1   +    6r)  (150^^02    -    3a3    -    16a^^) 

8tt 

8a-  3a,  3a-,  „    3a,         ^ 

-    p(15a,    ^—  +    15a^    -^  -3    -^   -    48a,       -^)  ]  T  ^ 
"^  13p  23p  3p  1        3p 


4 

+  1/2  (^^)       [(  1  +  7r)  (54a, '^  +  12a^^  -  75a,  ^a^  +  18a,  a,  -3a.) 
c  1  2  12  13  4 


3    3ot,  3a2  3ot-|  2    ^"^2 

-    p  (216a,       -TT —  +    24a-    -x —  -    150a,  a-    -x —  -  75a,       -^-— 
^              13p                  23p                     12     3p  13p 

3a^  3a,  3a.         „ 

+    18a,    -^  +    18a^   -x-f  -    3   -5-^)  ]  t"  }  (2.14) 

Idp  Jdp  dp 
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where  r  =  —  tt—  is  the  Gruniessen  constant.   This  expression 
c  dp  '^ 

is  the  first  four  terms  of  the  MR  formula.   The  standard  form 

2      4 
for  most  theories  is  e (p)  =  cp{l  -  yp   -  5p  ),  i.e., 

a-,  =  a..  =  0,  a-  ~  -  Y  and  a.  =  -6.   For  this  case  Eq.  (2.13) 

becomes 

_  _  TTV(27rk)^  ,4    10  ^,2jrk  2  6 
^  180h  c 

4 
+  1/2(12y^  +  36)  (^^)  T^]  (2.15) 

and  Eq.  (2.14)  reduces  to 

Pnh  -    --^4  ^(1  --  ^HT^-  +  i^  (^^^)^(l  +  SDy  -pfjlT^ 

+  3/2(2l]S)^(l  +  7r)(4Y2  +  6)  -  p  (8y|1  +  ||)  ]  T^} 


(2.16) 


As  a  check  on  the  validity  of  these  calculations  for 
pressure,  E,  C   and  S  were  calculated  from  the  free  energy — 
Eqns.  (2.9)  and  (2.15).   These  expressions  for  E,  C^  and  S 
were  compared  v;ith  other  calculations  of  those  quantities  by 
Zasada  and  Pathria  (1972)  ,  Roach  et   al .     (1972e)  and  Donnelly 
(1967)  and  were  found  to  be  in  exact  agreement.   These  calcu- 
lations led  to  the  discovery  that  Eq.  (2)  in  PVJH  (p.  1260)  is 
incorrect.   The  correct  form  is 

.  ViS^^   +  ^(21^)  V  +  7(4y2  +  5)(2^)V...] 
15h-^c''        '  ^ 


34 

Numerical  Results 
The  first  term  in  Eq.  (2.16),  tt  (27Tk)  ^  (1  +  3r)  T^^/lSOh-^c^  , 
common  to  all  theories,  will  be  referred  to  as  the  phonon 
pressure.   Similarly,  Eq.  (2.12)  gives  the  roton  pressure. 
To  compute  the  magnitude  of  the  pressure  due  to  the  disper- 
sion term  in  the  phonon  spectrum,  one  needs  the  values  of  y 
and  3y/8p.   Table  I  gives  the  values  of  y  and  9y/9p  for  three 
phonon  spectra  derived  from  experimental  results.   The  values 
of  y  and  3y/8p  for  most  of  the  theories  of  super fluid  helium 
fall  within  the  range  of  values  of  these  three  spectra.   The 
value  of  9y/9p  is  calculated  by  assuming  a  linear  variation 

of  y  with  p,  i.e.,  3y/3p  =  (y^^  ^^^  -  Ysvp)/(p24  atm.  "  ^SVP^ 

Because  this  thesis  concerns  the  measurement  of  pres- 
sure, numerical  estimates  of  the  size  of  the  various  contri- 
butions to  the  pressure  are  necessary.   Table  II  gives  numeri- 
cal values  of  the  phonon  pressure,  roton  pressure,  and  the 
pressure  due  to  dispersion  for  each  of  the  spectra  in  Table  I. 
These  values  are  given  at  saturated  vapor  pressure  and  at  24 
atm. — the  two  extremes  of  pressure  accessible  in  the  liquid. 
Values  for  both  1  K  and  0.5  K  are  given  since  the  upper  limit 
of  validity  for  the  approximation  that  there  are  only  phonons 
and  rotons  present  is  1  K,  and  since  the  effects  of  disper- 
sion are  clearest  near  0.5  K.   For  all  of  the  spectra  except 
MR  the  dispersion  pressure  varies  as  T  .   For  MR  the  disper- 
sion pressure  varies  as  T  ,  which  should  make  MR  readily  dis- 
tinguishable from  the  other  spectra,  although  the  absence  of 
any  comment  by  MR  on  the  variation  of  their  spectrum  with 
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Table  I.   Values  of  the  Dispersion  Parameter,  y,    and  Its 
Density  Derivative  in  c.g.s.  Units,  Based  on 
Experimental  Data 


Method 

SVP 

24  Atmospheres 

8Y 

9p 

Ultrasonic 

8  X 

,^37  a 

6  X  10^7  ^'^ 

-  7.3  X  lO^^ 

Neutron 

0  ^ 

6  X  10^7  d 

2.19  X  10^^ 

Specific  Heat 

-  4. 

1  X  1037  c 

19.6  X  10-^''  ^ 

8.65  X  10^^ 

^Eckstein  and  Varga  (1968) 

tiMills  et   at.     (1974),  Svensson,  Woods  and  Martel  (1972) 

^Cowley  and  Woods  (1971) 

^Svensson,  V7oods  and  Martol  (1972) 

^^Phillips,  VJaterfield  and  Hoffer  (1970) 
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Table  II.   Contributions  of  Different  Terms  to  the  Pressure, 
Expressed  in  Atmospheres,  at  Saturated  Vapor 
Pressure  and  at  2  4  Atmospheres 


Term  0.5  K  1.0  K 

1.47  X  10~^  2.35  X  10"^ 

3.33  X  lO"  5.34  x  10 

-1.18  X  10"^  -9.29  X  10~^ 

-3.29  X  10"^  -5.33  x  lO"'^ 

1.47  X  10"^  1.42  X  10"^, 

3.04  X  10  -4.79  X  10  ^ 


9.71  X  10  I  6.21  X  10_J 

9.74  X  10"^  6.23  x  10  ^ 


Phonon 
SVP 
24  atm. 

Roton 

SVP 

24  atm. 

Phonon  plus 
SVP 
2  4  atm. 

Roton 

Dispersion 

Ultrasonic 

SVP 

24  atm. 

Neutron 
SVP 
24  atm. 

Specific 
SVP 

24  a 

:  Heat 
itm. 

-2.27  X  10 J^  -1.48  X  10^ 

3.07  X  10"^  1.96  X  10  ^ 


-1.37  X  lO"'^  -8.78  X  10"^ 

7.7  X  10"'^  4.94  X  10"^ 
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density  makes  it  impossible  to  calculate  a  value  for  the  dis- 
persion pressure  for  their  phonon  spectrum.   The  values  of 
pressure  in  Table  II  should  be  compared  with  P  ,  which  ranges 
from  approximately  0.1  atm.  for  an  all-liquid  sample  just 
above  saturated  vapor  pressure  to  approximately  24  atm.  for 
an  all-liquid  sample  just  below  the  melting  curve. 

o 

It  is  difficult  to  calculate  the  T   term  in  Eq.  (2.16) — 

2 

the  next  order  correction  term — since  it  depends  on  y  , 

o 

Y  3y/9P/  6/  and  96/8p.   To  estimate  the  size  of  the  T   term 

at  1  K,  Y  and  3y/3P  were  taken  from  PWH.   6  was  set  equal  to 

2 

4y   and  8  6/8p  was  set  equal  to  By  3  6/3p,  so  that  the  contri- 

o 

bution  of  6,  and  86/8p  to  the  T  pressure  term  would  be  equal 

to  the  contribution  of  y  and  3y/9p.   This  gave  a  value  of 

76  78 

6  =  10    and  36/3p  =  -2  x  10    at  saturated  vapor  pressure. 

7  ft 

Only  Jackie  and  Kehr  (1974)  have  a  6  greater  than  10    at 
saturated  vapor  pressure.   Ho\-/ever,  the  large  positive  vari- 
ation of  their  6  with  p  would  tend  to  cancel  the  effect  of 

P 
the  large  5  term.   Despite  this  maximization  of  the  T   term, 

it  is  less  than  10%  of  the  T   term  at  1  K  and  decreases  more 

rapidly  v;ith  decreasing  temperature  than  the  T   term. 

It  is  also  necessary  to  estimate  the  error  caused  by  the 

inaccuracy  of  the  phonon-roton  approximation  in  the  region  of 

the  maximum  of  the  excitation  spectrum,  shown  in  Figs.  2  and 

3.   If  one  introduces  another  type  of  excitation,  called 

raaxons,  characterized  by  an  equation  of  the  form 

.  2 

,   (P  -  Po  ) 

e(p)  =  A r-^— 

2y 


38 


over  a  restricted  range,  the  properties  of  the  superfluid  can 
be  analyzed  in  terms  of  three  excitations  instead  of  two. 
Unfortunately,  because  of  the  necessity  of  restricting  the 
range  of  integration  on  each  of  the  excitations,  the  inte- 
grals in  Eqns.  (2.7)  and  (2.8)  can  no  longer  be  done  in 
closed  form.   They  must  be  done  numerically.   Since  the  spec- 
trum for  the  maxons  is  similar  to  the  roton  spectrum,  the 

equation  for  pressure  will  be  similar  to  Eq.  (2.10).   Thus 

I 

the  maxon  pressure  will  vary  as  exp(-AAT).  Substituting  a 

reasonable  A   in  Eq.  (2.11)  gives  a  pressure  equal  to  approxi- 
mately 5%  of  the  roton  pressure  for  low  pressures.   Since  the 
roton  pressure  is  considerably  smaller  than  the  phonon  pres- 
sure, the  effect  of  the  maxons  is  negligible  for  temperatures 
below  1  K.   As  shown  by  Fig.  3,  at  higher  pressures  the  situ- 
ation is  different,  the  height  of  the  maxons  changes  only 
slightly,  but  the  width  approximately  triples.   Because  A 
increases,  the  pressure  due  to  the  maxons  will  be  positive. 
A  reasonable  estimate  is  that  the  contribution  of  the  maxons 
will  be  approximately  10%  of  the  roton  pressure  at  higher 
density.   Because  the  roton  pressure  is  equal  to  or  larger 
than  the  phonon  pressure  at  these  densities,  this  will  intro- 
duce a  serious  error.   From  Fig.  3  it  would  appear  that  the 
Mills  spectrum,  derived  from  thermodynaniic  measurements,  is 
shifted  to  lov;er  p,  higher  A,  and  broader  y  in  order  to  aver- 
age in  the  effects  of  the  maxons.   Thus  the  difference  between 
the  thermodynamic  values  for  y,  p^  and  A  and  the  neutron 
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values  at  high  densities  may  be  due  to  the  inadequacy  of  the 
phonon-roton  approximation  in  the  region  of  the  maximum  of 
the  spectrum.   If  the  maxons  produce  a  pressure  large  enough 
to  interfere  v/ith  measurements  at  higher  densities,  one  could 
use  the  Mills  roton  spectrum  to  alleviate  the  problem. 

From  Table  II  it  is  evident  that  one  needs  a  pressure 

-5 
resolution  of  the  order  of  10    atm.  at  saturated  vapor  pres- 
sure to  clearly  see  the  dispersion  pressure.   At  24  atm.  a 
resolution  of  the  order  of  10~   atm.  is  required,  which  im- 
plies a  substantially  higher  instrumental  resolution.   If  such 
resolution  could  be  obtained,  pressure  measurements  would 
differentiate  between  the  different  predictions  of  the  dis- 
persion parameters. 


CHAPTER  III 
EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


The  Cryogenics 
The  cryostat  is  similar  to  the  one  described  by  Walsh 
(1963) .   It  has  been  substantially  modified  from  the  one 

described  by   Heberlein  (1969) .   The  main  features  of  the 

4 
cryostat  are  a  helium  bath,  a  vacuum  space,  a   He  refrigera- 
tor, a   He  refrigerator  and  a  sample  chamber.   Fig.  4  shows 
the  cryostat  schematically. 

The  vacuum  space  is  enclosed  by  a  copper  cylindrical 
container  which  is  attached  to  a  flange  on  the  cryostat.   The 

vacuum-tight  seal  is  made  with  0.075  cm  diameter  pure  indium 

4 
wire.   The  vacuum  space  surrounds  the   He  refrigerator,  the 

3 
He  refrigerator  and  the  sample  chamber.   When  the  apparatus 

is  first  being  cooled  down,  by  filling  the  heliiim  bath,  the 
vacuum  space  is  filled  with  helium  gas  for  purposes  of  ther- 
mal contact.   It  is  then  evacuated  by  pumping  for  several 
hours  with  a  Consolidated  Vacuum  Corporation  PMCS-2C  oil  dif- 
fusion pump  backed  by  a  Vfelch  model  1400  pump.   To  prevent 
contamination  there  is  a  nitrogen  trap  in  the  pumping  line. 

The   He  refrigerator  has  an  internal  volume  of  approxi- 

3 
mately  250  cm  .   It  is  filled  from  the  bath  via  a  modified 

Hoke  valve  operated  by  a  long  shaft  extending  through  the  top 

flange  of  the  cryostat.   This  refrigerator  can  be  maintained 
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Fig.  4.   Schematic  drawing  of  the  cryostat. 

A  Electrical  f eed-throughs 

B  Pumping  line  for  vacuum  space 

C  Diffusion  pump  for   He  refrigerator 

D  Mechanical  pump  for   He  refrigerator 

E  Trap  for   He  return  line 

F  Vacuum  space  flange 

G  He  refrigerator 

H  He  refrigerator 

I  Low  temperature  valve 

J  Radiation  shield 

K  Sample  chamber 

L  Vacuum  can 

M  Pumping  line  for   He  refrigerator 
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at  a  temperature  of  approximately  1.2  K  for  over  24  hours  by 
pumping  on  the  enclosed  helium  with  a  Kinney  model  KC-46 

pump . 

4 
A  radiation  shield,  attached  to  the  underside  of  the   He 

3 

refrigerator,  surrounds  the   He  refrigerator  and  the  sample 

chamber.   The  shield  was  constructed  in  the  following  manner. 
A  mat  of  #44  bare  copper  wire  was  wound  on  a  drum.   A  coat  of 
General  Electric  7031  varnish  was  applied  to  hold  the  mat 
together.   The  mat  was  then  fitted  around  and  epoxied  onto  a 
thin  phenolic  plastic  form.   The  upper  ends  of  the  copper 
wires  were  stripped  of  viirnish  and  soldered  onto  a  copper 

ring  which  screwed  into  another  copper  ring  bolted  to  the 

4 
bottom  of  the   He  refrigerator.   This  arrangement  provided 

the  largest  possible  cross-sectional  area  for  experimentation 

4 
below  the   He  refrigerator. 

3 

The   He  refrigerator  v/as  designed  for  continuous  opcra- 

3 
tion.   The  refrigerator  holds  approximately  1  cm  of  liquid 

3  3 

He.   To  increase  thermal  contact  between  the   He  and  a  cop- 

3 

per  flange,  which  is  the  base  of  the   He  refrigerator, 

approximately  1  meter  of  thin  (0.025  cm)  copper  foil  was 

placed  in  the  refrigerator  and  hard  soldered  to  the  base. 

3 

The  returning  He  gas  is  forced  to  pass  through  a  trap  im- 
mersed in  liquid  nitrogen.  This  trap  consisted  of  several 
layers  of  copper  wire  mesh  foD.lowed  by  Linde  molecular  sieve 

#13  X.   To  provide  the  pressure  drop  needed  to  liquefy  the 

3 
returning   He  gas,  an  impedance  was  placed  in  the  return  line 

below  the   He  refrigerator.   The  impedance  was  made  by 
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inserting  approximately  18  cm  of  9  mil  copper  wire  into  a 
10-mil-i.d.  copper-nickel   capillary.   A  National  Research 
Corporation  model  B-2  oil  diffusion  pump  backed  by  a  Welch 
model  1402  pump  v/ith  an  oil  shaft  seal  for  closed  system 
operation  was  capable  of  reducing  the  temperature  of  the 

sample  chamber  to  about  0.45  K.   By  closing  a  value  on  the 

3 

He  return  line  and  operating  for  periods  up  to  two  hours  in 

the  "single  shot"  mode,  colder  temperatures — about  0.39  K 
for  the  sample  chamber — v/ere  achieved. 

The  sample  chamber  is  supported  by  a  15  cm  long,  pitch 

bonded  graphite  rod  approximately  1  cm  in  diameter.   To  pro- 

3 
vide  thermal  contact  between  the   He  refrigerator  and  the 

sample  chamber,  a  copper  rod  with  a  diameter  of  0.5  cm  was 

3 
screwed  into  the  base  of  the   He  refrigerator  and  a  length 

of  copper  braid  with  a  diameter  of  0.4  cm  was  hard  soldered 

into  a  copper  rod  which  screwed  into  the  sample  chamber.   The 

two  copper  pieces  were  connected  by  a  lead  heat  switch.   The 

heat  switch  was  0.30  cm  wide  by  1  cm  long  by  0.0  3  cm  thick. 

It  was  joined  to  the  copper  pieces  v/ith  a  thin  coat  of  soft 

solder.   Around  the  switch  was  placed  a  small,  1  cm  bore, 

superconducting  solenoid  capable  of  generating  a  field  of 

1000  gauss  with  a  current  of  3  amps.   The  solenoid  was 

powered  by  a  6  volt  battery  which  also  provided  the  current 

for  the  heater  on  the  persistent  switch. 

The  Sample  Chamber  and  Pressure  Measurement 
The  sample  chamber  shown  in  Fig.  5  was  a  modification  of 
the  capacj  tive  strain  gauge  described  by  Straty  and  Adams 
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Fig.  5.   Sample  chamber 


A  Body  of  the  chamber 

B  Diaphragm  support 

C  Guard  ring 

D  Moveable  capacitor  plate 

E  Fixed  capacitor  plate 

F  Diaphragm 

G  Sample 

H  Inlet  Capillary 

The  copper  brush  is  not  shown 
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3 
(SA,  1969) .   The  cell  had  a  volume  of  approximately  5  cm  and 

a  height  of  0.4  2  cm.   To  increase  thermal  contact  to  the  sam- 
ple, the  interior  of  the  cell  was  filled  with  fine  copper 

wires  using  the  technique  of  Kirk,  Castles  and  Adams  (1971). 

2 
These  wires  had  a  surface  area  of  approximately  2  50  cm  and 

filled  about  30%  of  the  volume  of  the  cell.   The  main  body  of 
the  cell  was  made  of  copper,  while  the  diaphragm  and  the 
capacitor  plates  v/ere  made  of  berylium  copper.   The  diaphragm 
was  0.2  cm  thick.   The  upper  capacitor  plate  was  epoxied  onto 
a  post  machined  into  the  center  of  the  diaphragm,  so  that  as 
the  diaphragm  flexed  the  upper  plate  moved.   The  lower  capa- 
citor plate  was  held  in  a  fixed  position  by  a  guard  ring 
bolted  to  the  piece  containing  the  diaphragm. 

Changes  in  the  sample  pressure  cause  a  deflection  of  the 
diaphragm,  which  produces  a  change  in  the  capacitance  by 
changing  the  spacing  betv;een  the  capacitor  plates.   Expanding 
Eq.  (7)  in  SA,  one  can  show  that 

P-P   c.A(C-C)+B(C-C)^ 
o         o  o 

+  D(C  -  C  )^  +  E(C  -  C^)^  +  ...,         (3.1) 

where  P  is  pressure  and  C  is  capacitance.   In  practice,  the 
coefficients  D  and  E  are  so  small  that  the  contribution  of 
those  terms  is  negligible  except  for  very  large  pressures. 
Because  of  the  very  small  spacing  of  the  plates,  0.0012  to 
0.0025  cm,  any  imperfections  in  the  plates  will  cause  short- 
ing at  lov7  pressure.   The  plates  were  reasonably  free  of 
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imperfections,  since  a  pressure  of  about  15  atm.  was  reached 
before  shorting.   In  order  to  reach  higher  pressure,  the 
plate  spacing  was  increased  to  0.0075  cm.   With  this  spacing, 
the  plates  remained  unshorted  to  a  pressure  of  over  30  atm. 

The  capacitance  was  measured  with  a  General  Radio  type 
1620  A  capacitance  bridge.   An  Ithaco  model  391  A  lock-in  was 
used  to  detect  the  balance  of  the  bridge.   A  General  Radio 
1321  A  audio  oscillator  provided  a  10  v,  500  IIz  signal  to 
drive  the  bridge.   With  this  arrangement,  changes  in  capaci- 
tance as  small  as  10    pF  could  easily  be  detected.   This 
corresponds  to  a  pressure  change  of  from  1  x  10    atm.  to 
5  X  10~   atm. ,  depending  on  the  pressure.   As  the  pressure 
increases,  the  sensitivity  increases  because  the  plates  are 
closer  together.   Thus  this  apparatus  could  easily  detect  the 
dispersion  pressure. 

The  capacitance  was  measured  by  a  three-terminal  tech- 
nique which  uses  a  separate  coaxial  line  for  each  plate  with 
neither  plate  grounded.   This  prevented  shifts  in  capacitance 
due  to  small  movements  of  the  leads.   To  reduce  any  effect  of 
changing  temperature  on  the  bridge,  it  was  placed  in  a  large, 
insulated  box.   The  temperature  of  the  laboratory  was  regu- 
lated so  that  the  largest  temperature  drift  over  the  course 
of  a  run  was  1°  F.   The  measured  drift  rate  of  the  bridge 
never  exceeded  2  x  10~   pF  for  8  hours. 

The  calibration  of  the  strain  gauge  consisted  of  observ- 
ing the  capacitance  of  the  gauge  versus  the  pressure  read  on 
a  gauge  in  the  external  system.   The  external  gauge  was  a 
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Heise  model  #7770--a  dial  type  Bourdon  Gauge.   Calibration  of 
the  strain  gauge  extended  over  a  wide  range  of  pressures — 
from  0  to  10  or  14  atra.   The  calibration  points,  consisting 
of  a  capacitance  and  a  corresponding  pressure,  were  used  to 
find  the  coefficients  of  Eq.  (3.1)  by  means  of  a  least 
squares  program.   For  this  calibration  P  was  0  and  C  vvas 
the  capacitance  of  the  cell  at  4  K  when  evacuated.   The  re- 
sulting fit  had  an  RMS  deviation  of  less  than  0.01  atm.   The 
absolute  accuracy  of  the  calibration  measurements  was  no  bet- 
ter than  0.05  atm.,  the  absolute  accuracy  of  the  gauge.   How- 
ever, the  relative  accuracy,  i.e.,  the  accuracy  in  the  measure- 
ment of  small  changes  in  pressure  at  a  given  pressure,  was 
much  higher  being  limited  by  the  resolution  of  the  strain 
gauge  and  the  accuracy  of  the  derivative  dP/dC.   The  values 
of  dP/dC  must  be  fairly  accurate  to  obtain  a  good  fit  over  a 
wide  range.   The  thermodynamic  pressure  is  the  pressure  at 
constant  volume.   Obviously,  the  volunie  of  the  cell  is  not 
strictly  constant.   SA,  in  Eq.  (8) ,  show  that 

(9P)   =  dP       _1_  dV  ,   ^ 

^dT\         dT  ^^         k^V  dP^  ^^'     ' 

where  k_  is  the  isothermal  compressibility  of  the  sample. 
For  this  apparatus  the  correction  term  is  less  than  0.5%. 
Furthermore,  Boghosian  and  Meyer  (1966)  showed  that  k   is 
virtually  independent  of  temperature  below  1  K,  varying  by 
less  than  0.01%,  and  that  k   varies  very  slowly  with  pres- 
sure.   The  geometrical  factor  —  -^   is  independent  of  tempera- 
ture and  pressure  below  1  K  for  small  pressure  changes.   As 
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shown  in  Table  II,  the  total  variation  in  pressure  between  0 

-2 

and  1  K  is  approximately  10    atm.   Such  pressure  changes 

have  negligible  effect  on  k„  and  —  ^.   Thus  the  correction 
for  the  change  in  volume  of  the  cell  is  constant  with  tempera- 
ture below  1  K.   The  analysis  of  the  data  will  be  concerned 
mainly  with  the  temperature  dependence  of  the  pressure,  so  a 
small--less  than  0.5% — constant  correction  should  not  affect 
the  analysis  of  the  data.   For  that  reason  the  correction 
will  be  ignored  in  subsequent  analysis. 

Since  the  thermodynamic  pressure  is  the  pressure  at  con- 
stant volume  and  constant  density,  it  is  necessary  to  isolate 
the  cell  from  the  external  pressure  system.   At  first  a  valve 
on  top  of  the  cryostat  was  used.   This  proved  unsatisfactory 
becavise  of  variations  in  the  pressure  in  the  filling  capil- 
lary.  These  variations  were  caused  by  the  normal  changes  in 
the  level  of  the  liquid  helium  bath  which  changed  thermal 
gradients  in  the  filling  capillary  even  though  it  was  inside 
an  evacuated  tube.   To  alleviate  this  problem  a  low  tempera- 
ture valve,  shown  in  Fig.  6,  based  on  the  design  of  Roach  et 
at.     (1972d)  was  built.   The  valve  was  mounted  immediately  be- 
low the  sample  chamber  so  that  the  valve  was  at  the  same 
temperature  as  the  sample  cell.   The  valve  seal  v.'as  made  by 
a  Teflon-tipped  brass  stem  sealing  against  a  brass  seat.   The 
end  of  the  brass  stem  was  threaded  and  then  machined  to  a 
conical  shape.   The  inside  of  the  Teflon  tip  was  machined  to 
the  same  geometry  and  screwed  onto  the  brass  stem.   In  this 
way,  the  Teflon  tip  was  fully  supported  by  metal.   The  valve 
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Fig.  6.   Schematic  diagram  of  low  temperature  valve, 

A  Valve  body 

B  Accuator  bellows 

C  Valve  bellows 

D  Accuator  piston 

E  Valve  stem 

F  Valve  seat 
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was  an  hydraulic  one,  operated  by  changing  the  gas  pressure  in 
a  capillary  leading  to  the  valve.   At  the  start  of  an  experi- 
ment, helium  gas  was  let  into  the  capillary  and  liquefied  by 
contact  with  the  bath  until  the  valve  actuator  and  the  lower 
section  of  the  capillary  were  full  of  liquid  helium.   There- 
after, the  valve  could  be  opened  and  closed  easily  by  changing 
the  gas  pressure  in  the  filling  line.   A  pressure  of  about 
6  atjn.  was  required  to  fully  seal  the  valve.   A  small  SA  type 
strain  gauge  monitored  the  pressure  in  the  actuator. 

Unfortunately,  the  valve  had  a  small  superfluid  leak. 
Experimentation  showed  that  in  a  given  run  the  superfluid  leak 
rate,  i.e.,  -g^,  was  constant  in  time  regardless  of  pressure 
differential  across  the  valve.   No  dependence  of  the  leak  rate 
on  temperature  was  observed  for  temperatures  well  belov/  the  X 
line.   The  leakage  rate  changed  with  each  opening  and  closing 
of  the  valve,  varying  from  10    pF/min  to  10    pF/min.   Since 
the  leak  rate  was  constant  in  time,  a  simple  correction  fac- 
tor — the  e3apsed  time  multiplied  by  -^ — could  be  added  to 
each  capacitance  reading,  thus  removing  the  effects  of  the 
leak  on  the  capacitance  readings.   Upon  disassembling  the 
valve,  the  problem  appeared  to  be  small  flecks  of  solder  on 
the  valve  stem.   Two  subsequent  attempts  to  clean  the  valve 
to  prevent  leakage  failed,  resulting  in  much  larger  leaks. 

Temperature  Measurement  and  Regulation 
The  temperature  of  tlie  sample  chamber  was  measured  by  a 
calibrated  Cryocal  Germanium  therm.ometer  whose  electrical 
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resistance  varied  from  1000  ohms  at  the  X  line  to  400,000 
ohms  at  0.3  K.   The  thermometer  was  tightly  mounted  in  a 
copper  block  bolted  to  the  sample  chamber.   A  light  coat  of 
Apezion  'N'  grease  was  used  to  provide  thermal  contact.   The 
leads  to  the  thermometer  were  Midohm  wire  whose  resistance 

is  independent  of  temperature.   The  leads  were  thermally 

4  3 

anchored  to  the   He  refrigerator,  the   He  refrigerator  and 

the  sample  chamber  by  epoxying  them  to  copper  posts  bolted  to 
these  objects.   The  resistance  of  the  thermometer  was  mea- 
sured on  an  AC  resistance  bridge  described  by  Castles  (1973) . 
A  Princeton  Applied  Research  model  HR  8  lock-in  served  as 
both  the  source  of  the  20  Hz  driving  signal  and  as  the  de- 
tector of  the  balance  of  the  bridge.   A  G.R.   standard  decade 
resistor  v/ith  a  minimum  resolution  of  0.1  ohms  was  used  as 
the  known  resistance  in  the  bridge.   Using  this  arrangement 
it  was  easy  to  resolve  0.02%  changes  in  the  resistance.   This 
corresponds  to  a  temperature  resolution  of  better  than  0 . 1  mK 
over  the  entire  range  of  this  experiment. 

This  thermometer  had  been  calibrated  against  the  vapor 

3       4 
pressure  of   He  and   He  betv/een  0.4  K  and  4.2  K  by  Philp 

(1969).   He  found  that  a  form 

N       _  CIO  i-1 
T  =  I      A.  (R  •^■^^)  (3.3) 

i=l   ^ 

where  T  is  temperature  and  R  is  resistance  with  N  =  5  fitted 

the  temperature  "with  an  RI>1S  deviation  of  less  than  1  mK 

above  0.6  K."  (p.  39)   A  further  calibration  of  the  thermometer 
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against  the  susceptibility  of  CMN,  extending  to  0.29  K  was 
done.   These  results  were  fitted  by  the  same  form  using  the 
program  RESFIT  written  by  Philp.   These  fits  gave  a  set  of 
coefficients.  A.,  which  could  be  substituted  in  Eq,  (3.3)  to 
give  the  temperature  for  any  given  resistance.   This  yielded 
a  fit  with  an  RMS  deviation  of  less  than  2  mK  over  the  whole 
range,  with  a  slight  increase  in  the  deviations  near  0.3  K. 

The  out-of-balance  signal  of  the  lock-in  was  used  to 
modulate  the  power  supplied  to  a  heater  on  the  sample  cham- 
ber.  The  heater  v.'as  a  2000  ohm  metal  film  resistor  attached 
to  the  copper  braid  from  the  heat  switch,  just  above  the  sam- 
ple chamber.   A  metal  film  resistor  was  used  because  its 
resistance  was  almost  independent  of  temperature.   By  suit- 
ably adjusting  the  bias  current  and  the  magnitude  of  the  out- 
of-balance  signal  from  the  lock-in  with  an  external  circuit 
built  for  that  purpose,  the  temperature  of  the  sample  chamber 
could  be  stabilized  to  within  the  resolution  of  the  bridge  at 
any  temperature  betv/een  0.4  K  and  1.2  K.   Above  1.2  K  it  was 
necessary  to  reduce  the  cooling  power  of  the  refrigerator  be- 
cause the  heater  had  insufficient  power  to  stabilize  the 

3 

temperature.   This  could  be  done  by  pumping  on  the   He  v;ith 

3 
only  the  mechanical  pump,  or  by  turning  off  the   He  refrigera- 

3        4 
tor  entirely,  or  by  turning  off  both  the   He  and   He  refrig- 
erators, depending  on  the  temperatures  required.   For  the 
very  lowest  temperatures,  all  heating  and  regulation  were 
turned  off  and  the  system  was  allowed  to  slowly  drift  to 
colder  temperatures.   Thus  regulated,  the  temperature  was 
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quite  stable.   The  largest  observed  drift  was  4  mK  over 
several  hours. 

The  Sample 
The  sample  was  commercial  grade  helium  from  Air  Products 
and  Chemicals  Incorporated  (Airco) .   It  was  passed  through  a 

nitrogen  trap  with  the  same  copper  wire,  molecular  sieve  com- 

3 
bination  as  the   He  return  line  trap  previously  described. 

The  gas  then  passed  through  a  trap  filled  with  Linde  #13x 

molecular  sieve  immersed  in  liquid  helium.   Testing  with  a 

3 
Veeco  MS  9  leak  detector,  modified  to  scan  both  the   He  and 

He  peaks,  revealed  that  the   He  concentration  was  much  less 

than  0.01%. 

Procedure  for  Taking  Data 
Data  were  taken  by  recording  the  pressure  at  a  fixed 
temperature  after  the  system  came  into  equilibrium.   The  sam- 
ple v/as  then  warmed  or  cooled  to  a  new  temperature  and  ano- 
ther reading  was  taken.   In  most  runs,  the  apparatus  was 
first  allov/ed  to  cool  overnight.   Starting  from  this  tempera- 
ture, the  sample  was  systematically  warm.ed  in  approximately 
25  mK  steps  to  some  predetermined  temperature  greater  than 
1  K.   Then  the  sample  was  cooled  in  25  mK  steps  back  to  the 
original  temperature.   To  assure  consistency  between  runs, 
the  same  set  of  temperatures  was  used  for  each  data  pass.   In 
the  early  runs  that  extended  to  the  X  line,  it  was  found  that 
after  warming  to  the  X  line  the  apparatus  could  be  cooled 
well  below  the  starting  temperature,  but  would  suddenly  and 
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irreversibly  warm  up  to  above  the  starting  temperature.   This 
was  traced  to  inadequate  thermal  grounding  of  the  fill  capil- 
lary and  valve  actuator  capillary.   After  correcting  this 
problem,  it  v/as  possible  to  repeatedly  cover  the  whole  P  ver- 
sus T  curve  both  warming  and  cooling.   The  leakage  rate 
through  the  valve  was  measured  by  holding  the  cell  at  a  con- 
stant temperature  for  about  an  hour  and  observing  the  chance 
in  capacitance.   Since  the  P  versus  T  curve  has  a  maximum 
near  1.1  K  at  saturated  vapor  pressure  and  near  0.7  at  24 
atra. ,  the  leakage  rate  was  measured  at  this  point  to  minimize 
the  effect  of  any  inaccuracies  in  the  thermal  regulation. 


CHAPTER  IV 
RESULTS  AND  CONCLUSIONS 


Data  Reduction 


The  raw  data  consist  of  capacitance  values  versus  resis- 
tance bridge  readings  which  can  be  converted  to  pressure  ver- 
sus temperature.   Before  converting  the  capacitance  values  to 
pressures,  two  corrections  to  the  capacitance  were  made.   As 
described  in  Chapter  III,  the  correction  for  the  leakage 
through  the  valve  was  done  by  adding  a  constant  times  the 
elapsed  time  to  each  capacitance  value.   The  correction  is  al- 
ways additive  because  the  fill  capillary  is  evacuated  to  re- 
duce the  heat  leak,  so  that  the  mass  flow  is  out  of  the  cell. 
By  filling  the  capillary  to  a  pressure  greater  than  that  of 
the  cell,  the  mass  flow  is  into  the  cell,  the  capacitance 
rises  and  the  correction  is  subtracted  from  each  value.   In 

the  worst  case,  the  total  leakage  during  a  run  changed  the 

4 

density  of  the  liquid  by  less  than  one  part  in  10  .   Thus,  the 

approximation  of  constant  density  remains  valid.   After  apply- 
ing this  correction,  the  values  of  capacitance  for  both  warm- 
ing and  cooling  generally  coincided.   However,  in  some  of  the 
runs  there  were  obvious  discontinuities  in  the  capacitance 
values,  e.g.,  in  one  set  of  data  a  group  of  points  taken 
while  cooling  the  sample  v/as  exactly  4  x  10~   pF  larger  than 
the  same  points  taken  while  warning  the  sample.   Philp  had 
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reported  "Physical  shocks  .  .  .  caused  the  capacitance  to 
shift  discontinuously"   (p.  56).  The  problem  was  eventually 
traced  to  excessively  tight  vibrational  coupling  of  the 
cryostat  to  the  main  support  frame.   Small  contacts  with  the 
support  frame,  e.g.,  closing  the  -^He  return  line  valve  or 
bumping  an  exposed  corner  of  the  frame  would  repeatedly  cause 
discontinuous  shifts  in  capacitance,  although  moving  heavy 
objects  in  the  laboratory  or  banging  the  apparatus  violently 
other  than  on  the  main  frame  or  dewars  produced  no  shifts  in 
capacitance.   After  discovering  the  cause  of  the  problem, 
contact  with  the  frame  was  avoided.   No  further  discontinuous 
shifts  occurred. 

Corrections  for  the  shifts  were  applied  in  a  direction 
which  made  capacitance  versus  temperature  a  smooth  function. 
After  these  corrections  were  applied  to  the  capacitance 
values,  the  capacitances  taken  when  warming  the  sample  agreed 
very  closely  with  the  capacitances  taken  when  cooling  the 
sample.   In  all  cases  the  differences  between  warming  and 
cooling  were  within  three  of  the  smallest  resolved  units  of 
capacitance  over  the  entire  range  from  0.4  K  to  1.0  K  and 
were  frequently  zero.   This  agreement,  as  v/ell  as  the  sta- 
bility of  the  capacitance  at  a  given  temperature  indicated 
that  the  sample  was  in  thermal  equilibrium. 

The  corrected  capacitances  were  converted  to  pressure  by 
Eq.  (3.1),  using  the  coefficients  developed  in  the  calibra- 
tion procedure  described  in  Chapter  III.   The  resistances 
were  converted  to  temperatures  by  Eq.  (3.3)  using  the 
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coefficients  developed  in  the  calibration  procedure  described 
in  Chapter  III.   From  the  pressure  values,  the  density  was 
obtained  by  comparing  the  pressure  at  the  lowest  temperature, 
0.4  K,  with  the  pressure  versus  density  data  at  0.1  K  found 
in  Table  III  of  Abraham  et    al .     (Ab,  1970)  .   A  linear  interpo- 
lation was  used  between  their  data  points  which  introduced  a 
slight  error — no  more  than  0.0003  gra/cm   for  densities  between 
0.1451  and  0.1725  gm/cm-^ — because  the  pressure  versus  density 
function  is  not  strictly  linear.   To  within  the  accuracy  of 
the  data  of  Ab  (1970)  0.01  atm. ,  the  measured  pressure  is 
constant  below  about  0.8  K  so  that  no  significant  error  is 
introduced  by  using  the  pressure  value  at  0.4  K.   To  check 
this  determination  of  density,  measurements  were  extended  to 
the  X  line  and  the  pressure  at  the  \    line  was  compared  with 
the  measurements  of  pressure  versus  density  along  the  A  line 
by  Keesom  and  Keesom  (KK,  1933).   Good  agreement,  i  0.0005 
gm/cm  ,  was  found  between  densities  determined  by  comparison 
with  KK  and  those  determined  by  comparison  with  Ab  (1970)  , 
despite  differences  of  several  atmospheres  in  pressure  be- 
tween the  X  line  and  0.4  K  for  a  given  run.   Having  estab- 
lished this  correspondence  at  several  pressures,  subsequent 
runs  were  not  extended  to  the  X  line  because  of  experimental 
difficulties  in  doing  this.   For  this  reason,  the  values  of 
density  obtained  from  Ab  (19  70)  were  used. 

Knowing  the  density  of  the  sample  and  the  temperature  at 
each  point,  the  roton  pressure  can  be  calculated  from  Eq. 
(2.12).   Alternatively,  one  could  have  determined  the  roton 


parameters  from  the  data  and  compared  these  values  with 
values  determined  from  other  experiments.   This  was  not  done 
for  three  reasons.   (1)   Determining  the  roton  parameters 
from  the  data  would  introduce  six  new  parameters.  A,  3A/3p, 
\i,    d\i/dp,    p   and  ^Po/3p.   Four  of  these  parameters,   y,3y/3p, 
p   and  ^Po/9p,  are  difficult  to  measure;  hence  they  are 
poorly  known.   (2)   Such  an  analysis  would  increase  the  num- 
ber of  variables  to  10  for  the  simplest  spectrum.   This  is 
too  many  to  be  well  determined  from  just  the  20  data  points 
in  a  run.   (3)   This  would  have  required  non-linear,  least 
squares  fitting  routines  which  are  difficult  to  develop  and 
use,  and  can  lead  to  spurious  results. 

Since  P  =  P   +P   +?../  where  P  ,  includes  both  the 
o    r    ph         ph 

P  j^',  the  term  proportional  to  T  ,  and  the  dispersion  pres- 
sure, the  only  remaining  parameter  to  be  determined  is  P^, 
the  pressure  at  0  K.   Ideally,  this  would  be  determined  by 
lowering  the  temperature  until  the  pressure  reaches  a  con- 
stant value.   Since  this  requires  a  temperature  of  less  than 
0.25  K,  this  was  impossible  in  this  experiment.   Another, 
less  accurate,  possibility  is  to  extrapolate  the  lowest  tem- 
perature pressures  by  assuming  that  temperature  is  suf- 
ficiently low  that  the  only  contribution  to  the  pressure  is 

the  P  w'  which  varies  as  T^.   This  requires  a  substantial 
ph 

number  of  points  below  0.4  K.   Since  heat  leaks  limited  the 

experiment  to  T  greater  than  0.4  K,  such  a  method  was  useless, 

P  was  determined  from  a  least  squares  fit  of 
o 

P  -  P   =  P   +  AT^  +  BT^  +  CT^.   As  a  check  on  the  fit,  the 
r    o 
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sound  velocity  was  calculated  from  Ab  (1970) .   The  calculated 
value  agreed  closely,  r  1  to  3  m/sec  out  of  250  to  300  m/sec, 
with  the  experimental  values  given  in  Ab  (1970)  .   Other  in- 
ternal evidence,  discussed  later,  suggests  that  this  gave  a 

reasonable  value  for  P  .   The  reduced  data  are  a  set  of  values 

o 

ofPv,=P-(P   +P)  versus  temperature, 
ph         or 

Because  the  least  squares  routine  used  will  not  accept 
more  than  one  value  of  pressure  for  each  temperature,  the 
points  taken  cooling  and  warming  were  averaged  together  to 
produce  a  single  pressure  at  a  given  temperature.   The  co- 
efficients determined  by  the  least  squares  fits  were  converted 
to  meaningful  quantities  by  using  Eq.  (2.16)  with  the  values 
of  C,  r  and  p  from  Ab  (1970)  and  the  fundamental  constants  from 
Appendix  I  of  Donnelly  (1967) . 

Data  Analysis 
In  order  to  display  the  T^  and  T   terms  in  the  expres- 
sion for  P    as  given  by  Eq.  (2.16),  a  plot  of  P  i^/T   versus 

2 

T  was  made.   Fig.  7  shows  this  plot  for  the  experimental 

densities.   If  there  were  no  dispersion,  i.e.,  if  only  the 
phonon  pressure  were  present,  the  graph  v/ould  be  a  horizontal 
line.   On  this  type  of  plot  a  T   dependence  is  a  straight, 
non- horizontal  line.   The  lines  along  the  left  edge  of  the 
figure  indicate  the  position  of  this  horizontal  line  for  each 
run.   The  lines  are  labeled  with  the  density  of  the  sample  for 
that  run.   This  is  a  plot  of  the  "raw"  data,  i.e.,  every  ob- 
served point  in  this  temperature  range  has  been  reduced  and 
plotted.   At  higher  temperatures  the  difference  between  points 
is  smaller  than  the  size  of  the  symbol  used  to  plot  them.   The 


61 


0,022 


0.020 


0.03.8 


0.016 


si 
0, 
^        0.014 


0.012 


0.010 


0.14G4 


0.1474 


0.1496 


0.1519 


0.1549 


=iLJ576 


O.i 


«o 


4 


•  3 


T 


O 

o 


o  o  o 


o 
o 


n  a 
a  n 


V  f 


T     V     V 


•   «    « 


o    o   o 
o 

o 


O  O 


o  o 


■    a 


L 


0.2 


0.3 

m2 


0.4 


0.5 


0.6 


(k2) 


Fig.  7.   Pph/T'^  versus  T^.   Lines  indicate  position  of  line  of 
zero  dispersion  for  density  given  on  the  line.  Solid  tri- 
angle 0.1464,  solid  circle  0.1474,  open  circle  0.1496, 
solid  square  0.1519,  open  triangle  0.1549,  open  square 
0.1576. 
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scatter  at  the  lower  temperatures  is  largely  due  to  instru- 
mental resolution,  magnified  by  the  T~^  weighting  factor. 
The  error  bars  represent  instrumental  resolution.   The  smooth- 
ness with  which  the  plots  approach  the  zero  of  temperature  is 

indicative  of  the  correctness  of  the  method  of  choosing  P  . 

o 

Experimentation  showed  that  small  changes  in  P   produced  pro- 
nounced curvature  in  the  low  temperature  end  of  the  plot — 
upward  when  P  was  too  low,  downward  when  P  was  too  high. 
In  the  low  temperature  limit,  the  dispersion  pressure  becomes 

negligible  so  that  P    is  equal  to  P  ,  .   Thus  these  plots 

ph  ph 

must  join  the  horizontal  lines  smoothly  in  the  low  tempera- 
ture limit.   The  two  highest  densities  are  smooth  in  the  low 
temperature  limit  hut  do  not  extrapolate  to  the  phonon  pres- 
sure.  The  reason  for  this  is  unknown. 

The  most  noticeable  and  unexpected  feature  of  these  plots 
is  the  sharp  upward  curvature  of  each  plot  as  the  temperature 
increases.   Because  the  curvature  is  upward,  this  represents 
an  "excess  pressure."   Since  the  excess  pressure  contribution 
has  a  notice^lble  convex  curvature,  its  temperature  dependence 
is  stronger  than  T  .   For  higher  densities,  the  magnitude  of 
the  excess   pressure  is  larger  relative  to  P  ,  .   Also  at 
higher  densities,  the  excess  pressure  becomes  noticeable  at 
lower  temperatures.   The  curvature  of  the  plots  due  to  the 
excess  pressure  contribution  is  greater  at  higher  densities 

than  it  is  at  lower  densities. 

g 
To  test  whether  this  excess  pressure  arose  from  the  T 

term  in  Eq.  (2.16),  the  points  up  to  T  --  0.85  K  were  fitted 
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with  a  least  squares  routine  to  the  functional  form 
P  =  Pj^  +  AT   +  BT^  +  CT   +  P  .   Then,  systematically  the 
high  temperature  points  were  removed.   This  sharply  lowered 
the  RMS  deviations,  typically  by  a  factor  of  five  when 
changing  from  20  to  14  points,  and  substantially  changed,  by 
a  factor  of  two,  the  value  of  the  coefficients  B  and  C. 
Furthermore,  the  values  of  B  and  C  for  points  up  to  0.8  5  K 
were  well  outside  the  range  of  predicted  values.   In  contrast, 
the  values  for  B  and  C  for  points  up  to  0.7  K  were  in  reason- 
able accord  with  predicted  values.   The  magnitude  of  the  ex- 
cess pressure  contribution  substantially  exceeds  the  value 
estimated  for  the  T   term  in  Chapter  II.   Thus  it  was  con- 

g 

eluded  that  excess  pressure  v;as  not  due  to  the  T   term  in  Eq. 
(2.16)  . 

The  higher  order  terms  in  Eq.  (2.16)  and  in  the  expan- 
sions of  other  thermodynamic  functions  are  small  since  these 
expansions  are  rapidly  convergent  with  increasing  pov/ers  of 
T,  so  it  seems  unlikely  that  the  excess  pressure  arises  from 
omitted  higher  order  terms.   While  there  is  some  evidence 
that  the  shape  of  the  spectrum,  hence  y  and  5,  may  vary  with 
temperature,  a  variation  with  temperature  of  the  magnitude 
required  to  produce  excess  pressure  seen  here  would  totally 
destroy  the  Landau  approach  which  is  well  supported  by  other 
measurements.   The  roton  pressure  could  be  in  error,  but 
there  is  no  reason  to  believe  that    P   is  incorrect. 

The  likeliest  explanation  for  the  excess  pressure  is  the 
maxons,  i.e.,  contributions  from  the  maximum  of  the  excitation 
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spectrum  which  are  poorly  approximated  by  the  phonon-roton 
model.   As  was  argued  in  Chapter  II,  the  maxon  contribution 
to  the  pressure  would  be  positive,  increase  with  increasing 
density  and  would  vary  as  exp  (-l/T) .   This  is  consistent 

ACQ 

with  the  observations,  since  T   +  T   +  T  V70uld  be  a  poor 
approximation  to  exp  (-1/T) .   In  addition,  there  is  evidence 
that  a  similar  phenomenon  was  observed  in  the  specific  heat 
measurements.   In  their  analysis,  PWH  had  to  limit  themselves 
to  T  less  than  0.7  K  to  achieve  good  fits.   Furthermore,  the 
PWH   value  of  Y  at  20  atm.  is  a  factor  of  3  higher  than  y 
determined  by  neutron  scattering  and  ultrasonic  measurements, 
which  are  in  good  agreement  with  each  other.   PWH  let  the 
data  determine  the  roton  parameter,  A,  in  contrast  to  this 
work  which  calculates  the  roton  contribution  to  the  pressure. 
The  PWH  values  of  A  agree  much  more  closely  v.'ith  Mills  (1965) 
values  than  with  Donnelly's  (1972)  values.   As  noted  in  Chap- 
ter II  the  Mills  spectrum  disagrees  with  the  neutron  scatter- 
ing values  in  a  fashion  that  would  suggest  an  attempt  to 
average  over  the  maxons.   The  Mills  spectrum  is  based  on 
thermodynamic  measurements.   This  suggests  that  all  thermo- 
dynamic data  see  the  effect  of  maxons  but  that  these  contri- 
butions are  obscured  by  allowing  the  roton  parameters  to  vary, 
Even  with  this  variation  of  the  roton  parameters,  the  Mills 
spectrum  remains  a  poor  approximation  to  the  maxons.   Such  an 
error  in  approximation  might  explain  the  discrepancy  betv;een 
Y  of  PWH  at  20  atm.  and  y  from  other  measurements.   Neutron 
scattering  determines  y  by  fitting  the  spectrum  in  the  small 
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p  region.   This  fit  is  cut  off  below  the  maximum  of  the  spec- 
trum.  Thus  the  maxons  have  no  effect  on  the  determination  of 
Y-   Ultrasonic  measurements  are  usually  taken  at  0.1  K,  well 
below  the  temperatures  at  which  maxons  contribute.   Both  the 
specific  heat  and  these  pressure  measurements  extended  to 
0.9  K.   At  that  temperature,  the  calculation  of  the  thermo- 
dynamic properties  directly  from  the  excitation  spectrum  at 
saturated  vapor  pressure  by  Singh  (1968)  showed  that  there 
were  detectable  discrepancies  between  calculations  based  on 
the  phonon-roton  model  and  direct  calculations.   It  appears 
from  these  data  that  these  discrepancies  are  larger  and  be- 
come detectable  at  lower  temperatures  as  the  density  is  in- 
creased. 

Turning  to  the  lov;  temperature  portion  of  the  data,  it 
is  evident  that  all  of  the  curves  have  an  initial  downward 
slope.   This  slope  appears  to  decrease  with  increasing  density. 

The  coefficient  of  the  T   term  in  Eq.  (2.16)  is  proportional 

2 

to  1/c  ,  where  c  is  the  sound  velocity.   Since,  as  the  density 

increases,  c  increases,  some  decrease  in  slope  of  the  plots  in 

Fig.  7  would  be  expected.   To  calculate  the  value  of 

y'  ^  (Y  -  i^^rp  |-^)  from  Eq.  (2.16)  at  different  densities, 

J.+  ->  1   dp 

the  following  procedure  was  adopted.   Values  of  c  and  T,    taken 

from  Ab  (1970)  were  used  to  calculate  P  ,  '  and  the  constants 

fi       8 
in  the  coefficients  of  the  T   and  T   terms  in  Eq.  (2.16).   The 


quantitv  P-(P   +P   +P,')  was  least  square  fitted  to  both 
^  '  or    ph 

fi  I  Pi       8 

a  form  BT   and  a  form  B  T   +  CT  .   Data  points  were  systema- 
tically removed  from  the  high  temperature  end  until  the  effects 
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of  the  excess  pressure  were  negligible.   Two  criteria  were 
used  to  determine  the  absence  of  the  excess  pressure.   (1) 
The  quality  of  the  fit  from  the  data  for  the  form  BT  must 
not  depend  systematically  on  T.   If  the  data  used  in  the  fit 
extend  into  the  excess  pressure  region,  the  deviations  in- 
crease sharply  at  the  high  temperature  end.   (2)   The  value 
of  Y   derived  from  B  must  agree  reasonably  well  with  the  value 
derived  from  B  . 

The  results  of  this  procedure  are  shown  in  Table  III  for 
the  four  lowest  densities  for  both  the  T   fit  and  the  T   +  T 

fit.   The  three  higher  densities  fail  to  join  the  horizontal 

4 
AT   line  smoothly  at  T  =  0,  probably  due  to  the  increasing 

magnitude  of  the  excess  pressure,  so  that  fits  of  this  type 

are  useless.   The  actual  numbers  are  accurate  to  only  ±  10  to 

25%  becavase  of  the  relatively  small  number  of  points  fitted 

(10  to  15)  and  the  sizeable  changes  that  occur  when  the 

points  are  removed.   Nevertheless,  it  is  clear  that  y   is 

negative  and  decreases  in  magnitude  as  the  density  increases. 

The  numbers  for  6*  =  4y^  +  6-m-r  (8y  |^  +  P-)  from  Eq.  (2.16) 

IH-  /!        op     op 

are  about  equal  to  4y   -  \?nv    (^y  --n— )  /  the  contribution  of 

1+ 11  dp 

higher  powers  of  y-   There  is  very  little  improvement  in  the 

o 
RMS  deviations  when  the  T   term  is  included.   These  values  of 

Y   should  be  compared  v/ith  values  of  y   based  on  the  spectra 

37 
m  Table  I.   At  saturated  vapor  pressure,  these  are  8.5x10 

37 
for  the  ultrasonic  case,  -2.5  x  10    for  the  neutron  case, 

38  ' 

and  -1.26  x  10    for  the  specific  heat  case.   A  plot  of  y 


versus  density  is  shown  in  Fig.  8 
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I       I 

Table  III.   The  Values  of  y   and  5   at  Experimental  Densities, 


Density  ,a               ,b               .b 

in  gr/cm-^      Y Y o 

0.1464  -1.06  X  10^^  -1.75  x  10^^  8.75  x  lo''^ 

0.1472  -1-19  X  10^^  -1.56  x  10^^  5.10  x  10'^^ 

0.1496  -9.95  X  10^"^  -1.95  x  10^^  9-84  x  lo"^^ 

0.1519  -5.80  X  10^"^  -6.64  x  10^'  1.94  x  lo''^ 


^derived  from  the  T   fit 
^derived  from  the  T^  +  T^  fit 
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8.   y'  versus  density.   Open  triangle  this  work,  open 
square  PWH  assuming  Y  of  CW  at  24  atm.  ,  open  circle  PIVH, 
closed  triangle  ZP  assuming  y  of  CW  at  24  atm. ,  closed 
square  ZP,  dashed  line  is  based  on  normal  dispersion, 
dashed  and  dotted  line  is  based  on  CW. 
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To  illustrate  the  agreement  of  the  pressure  data  with 
other  data.  Fig.  9  shows  a  comparison  of  the  best  data  with 
a  calculated  plot  for  the  same  density,  assuming  y   has  the 
PWH  value  as  corrected  by  Zasada  and  Pathria  (1974)  and  cal- 

culating  |^  =  (734  ^tm.  -^^SVP^/^P24  atm.  "  ^SVP^  '  ^^^^^^^9 

Yo„   ..    has  the  CW  value  and  Jc-^tt,   has  the  P^'TH  value.   The 
'24  atm.  SVP 

data  clearly  follow  the  trend  of  the  calculation  up  to  the 
temperatures  where  the  excess  pressure  becomes  detectable. 
Further  jntercomparison  of  theories  is  shown  in  Fig.  9  by 
plotting  y'  derived  from  the  BT   fit  versus  density.   For 
comparison,  the  values  of  y   versus  density  for  the  specific 
heat  data  are  shown.   These  are  calculated  using  values  of  y 
of  PWH  and  those  of  PWH  as  corrected  by  Zasada  and  Pathria. 
The  value  of  8y/9p  is  calculated  by  assuming  a  linear  fit  be- 
tween the  lowest  and  highest  densities  of  PWH.   Also  shown 
are  the  values  of  y   based  on  the  ultrasonic  and  neutron  spec- 
tra given  in  Table  I.   It  is  evident  that  there  is  good  agree- 
ment between  values  of  y   derived  from  pressure  measurements 
and  those  derived  from  specific  heat  measurements;  and  there 
is  poor  agreement  with  the  neutron  and  ultrasonic  values. 
The  ultimate  goal  of  these  measurements  is  to  derive 
values  of  Y  as  a  function  of  density.   To  do  this  one  takes 
the  values  of  y   as  a  function  of  density  and  calculates  y    in 
a  self-consistent  way.   This  is  done  by  choosing  an  initial 
value  of  Y  and  dy/dp   which  yield  the  correct  y      at  that  den- 
sity.  These  are  used  to  calculate  y    at  another  density,  which 
in  turn  determines  3y/9p  at  that  density  from  y    .   These 
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values  are  used  to  calculate  y  and  oy/dp    at  the  next  density. 
The  process  is  repeated  until  a  consistent  set  of  y  and  Sy/Sp 
are  calculated,  i.e.,  dy/dp    and  y  are  smooth  functions  of  den- 
sity.  For  the  present  data,  this  is  a  very  imprecise  calcula- 

I  I 

tion  because  there  are  too  few  values  of  y  ,  the  values  of  y 

are  too  uncertain  and  the  range  of  densities  covered  is  too 
restricted. 

If  the  excess  pressure  contributions  were  removed,  there 

would  be  a  much  greater  range  of  P  versus  T  from  which  to 

I  I 

extract  y  .   This  would  greatly  improve  the  accuracy  of  y  . 

Also,  y   could  have  been  determined  over  a  wider  range  of  den- 
sities.  Nevertheless,  the  fact  that  y   is  negative  and  in- 
creases with  increasing  density  leads  to  the  conclusion  that 
at  low  pressure  y  is  negative  and  that  9y/9p  is  positive. 

The  possibility  of  a  linear  dispersion  term  was  conjec- 
tured by  MR.   To  check  this  possibility,  a  plot  for  the  best 

data  of  P  , /T   versus  T  was  made.   The  results  are  shovm  in 
ph^ 

5 
Fig.  10.   If  P  ,  varies  at  T   at  low  temperatures,  the  plot 

would  be  a  straight  line.   The  curvature  of  the  plot  at  low 

temperatures  is  obvious.   Furthermore,  the  curve  does  not 

extrapolate  properly  to  T  =  0.   This  contrasts  with  the  linear 

appearance  of  the  plot  of  P  ,  /T   versus  T   in  Fig.  8.   Thus 

the  dispersion  pressure  appears  to  vary  as  T  ,  rather  than  T  , 

at  low  temperatures.   From  Eq.  (2.14),  this  implies  that 

[  (l+4r)  (-aj-p-^]  =  0.   Since  all  of  the  plots  in  Fig.  7 

X      dp 

appear  to  be  similar  to  the  best  data  at  low  temperatures, 
one  can  conclude  that  3"i/8p  =  0  and  thus  a-j_  =  0.   This  con- 
tradicts the  MR  form  of  the  dispersion  curve. 
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Conclusions 

Three  conclusions  can  be  drawn  from  this  work.   (1)   The 
MR  form  of  the  excitation  spectrum  is  probably  incorrect. 
(2)   There  is  anomalous  dispersion  in  the  excitation  spectrum 
of  superfluid  helium  at  low  densities.   The  magnitude  of  the 
dispersion  parameter  and  density  derivative  have  values  simi- 
lar to  those  derived  from  the  specific  heat  measurements  of 
PV7H.   (3)   There  are  contributions  to  the  pressure  other  than 
phonons  and  rotons,  probably  maxons.   These  contributions  in- 
crease in  magnitude  as  the  density  is  increased.   These  con- 
tributions must  be  removed  before  accurate  conclusions  can  be 
made  about  the  values  of  y  and  Sy/Sp  as  a  function  of  density, 
particularly  at  higher  densities. 

More  work  needs  to  be  done  in  several  areas  before  pres- 
sure measurements  can  be  used  to  discriminate  among  the  dif- 
ferent theoretical  forms  for  the  anomalous  dispersion.   At 
present  the  uncertainties  in  this  experiment  are  too  great  to 
permit  this  to  be  done  with  any  confidence.   The  data  need  to 
be  improved  with  more  runs  of  the  quality  of  the  best  one. 
The  data  should  be  extended  to  lower  temperatures  if  possible 
to  permit  an  unequivocal  assignment  of  P   based  on  the  data. 
There  is  a  need  for  a  greater  number  of  data  points.   The 
25  mK  spacing  was  used  so  that  the  change  in  capacitance  for 
each  change  in  temperature  was  clearly  resolvable.   However, 
an  additional  set  of  data  could  be  taken  25  mK  apart  with 
each  point  being  the  middle  of  one  of  the  present  intervals. 
This  would  double  the  number  of  data  points  and  improve  the 
accuracy  of  the  fits. 
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A  method  of  reraoving  the  excess  pressure  should  be  de- 
vised so  that  the  full  range  of  temperature  can  be  used  in 
fitting  the  temperature  dependence  of  the  dispersion  term. 
This  will  require  an  extensive  neutron  scattering  study  of 
the  region  of  the  maximum  of  the  excitation  curve  as  a  func- 
tion of  density.   Such  work  is  in  progress  at  the  Brookhaven 
National  Laboratory.   Using  those  data,  one  could  calculate 
the  contribution  due  to  the  maxons.   Such  calculations  could 
be  used  to  reanalyze  both  these  pressure  data  and  the  speci- 
fic heat  data  of  PWH.   Such  a  reanalysis  would  improve  the 
accuracy  and  the  reliability  of  the  values  of  Y/  5  and  their 
density  derivatives. 

Studies  of  the  dispersion  curve  have  widespread  appli- 
cability, and  are  not  just  limited  to  studies  of  helium.   For 
example,  the  differences  among  the  theories  largely  arise 
from  differing  treatments  of  the  interactions  among  excita- 
tions.  Excitation  models  are  common  throughout  physics,  par- 
ticularly in  solid  state  physics  v/here  the  concept  of  phonons 
originated.   A  common  problem  in  all  of  these  models  is  the 
treatment  of  the  interactions  of  excitations .   In  a  solid  the 
problem  is  particularly  complicated  because  of  the  existence 
of  three  branches  of  the  excitation  curve.   Helium,  a  simple 
liquid,  is  a  good  substance  for  testing  different  theories. 
Thus,  studies  on  helium  could  yield  theories  applicable  to 
many  materials. 

Dispersion  curves  appear  not  only  in  solids  but  also  in 
liquids  and  amorphous  solids.   Goda  (1972)  has  noted  that 
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"longitudinal  phonon  dispersion  of  some  amorphous  solids  and 

simple  liquids  are  of  the  phonon-roton  type  as  observed  in 

4  ... 

liquid   He.   Such  behavxor  seems  to  exist  also  in  liquid 

metals"   (p.  1064).  Is  the  observed  anomalous  dispersion 
unique  to  helium,  or  is  it  a  common  feature  of  all  topologi- 
cally  disordered  substances?   Perhaps  precision  pressure 
measurements  of  liquid  neon  or  liquid  argon  could  answer 
that  question. 

The  universality  of  phonon-roton  type  dispersion  curves 
suggests  another  reason  for  making  detailed  studies  of  the 
dispersion  curve  of  superfluid  helium.   Khalatnikov  (1965) 
gave  a  simple  argument  that  claims  to  show  that  any  liquid 
with  a  roton  minimum  will  be  superfluid  at  T  =  0.   That  im- 
plied that  the  form  of  the  dispersion  curve  "caused"  super- 
fluidity.  The  Khalatnikov  argument  does  not  predict  any 
clearcut  transition  from  superfluid  to  normal  fluid  as  the 
temperature  is  raised  unless  the  spectrum  changes.   But 
helium  has  a  phonon-roton  spectrum  above  the  X  line,  and  is 
not  a  superfluid  at  that  point.   There  are  qualitative  dif- 
ferences in  the  excitation  spectrum  of  helium  above  and  be- 
lov7  the  X    line.   The  line  v/idth  is  much  broader  and  the 
temperature  dependence  of  quantities  such  as  the  roton  mini- 
mum are  different.   Takeno  and  Goda  (1972)  present  a  model 
which  successfully  predicts  the  spectrum  of  He  I  as  well  as 
other  simple  liquids  but  fails  for  He  II.   Accurate  knowledge 
of  the  excitation  spectrura  for  He  II  could  lead  to  a  better 
understanding  of  the  fundamental  nature  of  superfluidity. 


76 


This  work  has  shown  that  the  excitation  spectrum  of  super- 
fluid  helium  has  anomalous  dispersion  at  low  densities  and 
that  for  precision  measurements  the  contributions  of  the  maxons 
must  be  included  in  the  theory;  but  it  does  not  provide  the 
final  answer  as  to  which  theory  is  correct. 
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